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ABSTRACT 

We present near-infrared (NIR) color-magnitude diagrams (CMDs) for the resolved stellar popula- 
tions within 26 fields of 23 nearby galaxies (< 4Mpc), based on images in the i^llOPF and F160VF 
filters taken with Wide Field Camera 3 (WFC3) on the Hubble Space Telescope (HST). The CMDs are 
measured in regions spanning a wide range of star formation histories, including both old dormant and 
young star-forming populations. We match key NIR CMD features with their counterparts in more 
familiar optical CMDs, and identify the red core Helium burning (RHeB) sequence as a significant 
contributor to the NIR flux in stellar populations younger than a few 100 Myrs old. The strength of 
this feature suggests that the NIR mass-to-light ratio can vary significantly on short timescales in star 
forming systems. The NIR luminosity of star forming galaxies is therefore not necessarily proportional 
to the stellar mass. We note that these individual RHeB stars may also be misidentified as old stellar 
clusters in images of nearby galaxies. For older stellar populations, we discuss the CMD location of 
asymptotic giant branch (AGB) stars in the HST filter set, and explore the separation of AGB sub- 
populations using a combination of optical and NIR colors. We empirically calibrate the magnitude 
of the NIR tip of the red giant branch (TRGB) in FIQQW as a function of color, allowing future 
observations in this widely adopted filter set to be used for distance measurements. We also analyze 
the properties of the NIR RGB as a function of metallicity, showing a clear trend between NIR RGB 
color and metallicity. However, based on the current study, it appears unlikely that the slope of the 
NIR RGB can be used as an effective metallicity indicator in extragalactic systems with comparable 
data. Finally, we highlight issues with scattered light in the WFC3, which becomes significant for 
exposures taken close to a bright earth limb. 

Subject headings: galaxies: stellar content — galaxies: abundances — galaxies: distances and redshifts 
— galaxies: dwarf — galaxies: irregular — stars: AGB and post- AGB — stars: 
carbon — infrared: stars 



1. INTRODUCTION 



Near-infrared (NIR) observations have become increas- 
ingly important for studies of galaxies and their ev olution 
over cosmic time (e.g .,^onsclice et al. 2005; Dahlen et al. 



2005| [Saracco et al . 2006; Cirasuolo erar.]|2010[ ). I'his 



trend will no doubt continue during the coming decade, 
thanks to continued improvements in NIR detectors, the 
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growing maturity of adaptive optics, the installation of 
WFC3 on HST, and the upcoming launch of the James 
Webb Space Telescope (JWST). 

The importance of NIR observations for galaxy evolu- 
tion studies rests on the reduced sensitivity of the NIR 
mass-to-light ratio to dust and to the age or metallicity of 
the underlying stellar population, particularly compared 
to the optical or ultraviolet. The NIR luminosity of a 
galaxy is therefo re thought to be a robust indicator of 
stellar mass (e.g. 'T hronson k, Gree nho use | |1988 Brinch- 



mann & Ellis 2000;"Bundy et al. '5DD5). However, the 



ability to correctly interpret NIR observations relies on 
accurate stellar population modeling, which in turn re- 
quires accurate isochrones, stellar lifetimes, and spectra 
for the evolving stars which dominate the fiux at NIR 
wavelengths (see Conroy & Gunn (2010), for an estimate 
of current uncertainties) . 

Recently, there has been some concern that the mass- 
to-light ratio in the NIR is not nearly as stable as has 
been assumed. On the theoretical side, IMaraston et al.l 



( 2006 1 have suggested that the NIR fiux from asymptotic 



giant branch (AGB) stars may lead to drastically lower 
mass-to-light ratios when intermediate-age populations 
are present, as must be the case for in situ observations 
of young ga laxies at high redshift (although see Kriek| 
etal]p010| )). Lower NIR mass-to-light ratios could also 
help to explain why some elliptical galaxy progenitors at 
high redshift appear to have stellar densities higher than 
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seen in the local universe (van Dokkum et al][2008| ) 



The most accurate constraints on evolving stars' con 



ulations directly 


(e.g. Rejkuba et al. 


2006 


Gullieuszik 


et al. 2007 


2008 


Melbourne et al.|201L 


a most recently). 



solve large numbers of individual stars in nearby galaxies. 
There are well over a hundred galaxies within the Local 
Volume {D < A Mpc) that are close enough to be targeted 
for resolved stellar population studies. These galaxies 
have alr eady bee n observed extensively with HST in the 
optical dDaicanton et al.||2009[ [Gilbert et al.pOllj ), and 
their color-magnitude diagrams have revealed the galax- 
ies to have a wide range of metallicities and star forma- 
tion histories. These systems therefore form an ideal set 
for probing the properties of the same stellar populations 
in the NIR. 

To this end, we undertook an HST "snapshot" survey 
of nearby galaxies with existing high-quality optical data. 
We took advantage of the high throughput and resolv- 
ing power of the WFC3 IR channel to produce CMDs 
of individual fields in 23 galaxies, spanning a variety of 
stellar populations. In this paper, we present our sam- 
ple, and characterize the resulting color-magnitude di- 
agrams. These data will be used in subsequent papers 
to place quantitative constraints on the con tribution of 



AGB and R HeB stars to the total luminosity (Melbourne 
et al.|20ll|), on the lifetimes of AGB stars, and on the 



population of carbon stars. 

The outline of the paper is as follows. In Section [2] we 
present our sample selection, observations, reductions, 
and matched optical photometry. We also discuss a pre- 
viously uncharacterized transient scattered light feature 
in the Wide Field Camera 3 infrared channel (hereafter 
WFC3/IR). In Section [3] we discuss the origin of differ- 
ent features in the NIR CMDs and optical-NIR color- 
color diagrams, and highlight the importance of core red 
Helium-burning (RHeB) stars to the overall luminosity, 
for populations younger than ~ 0.5 Gyr. We discuss the 
properties of the tip of the red giant branch (TRGB) in 
the NIR (Section [4 1, and the AGB and red giant branch 
(RGB) luminositjAfunction in Section [sj Finally, we dis- 
cuss the metallicity dependence of the NIR RGB mor- 
phology and color on metallicity in Section [6j 



2. DATA 

2.1. Sample Selection 

To characterize NIR CMDs as a function of stellar age 
and metallicity, we require targets whose stellar popula- 
tions have been constrained using optical observations in 
tandem with well-calibrated stellar isochrones. We there- 
fore chose SNAP targets from among those nearby galax- 
ies with archival high-quality multi-color imaging (typ- 
ically from ACS, but in some cases from WFPC2). Of 
these, we selected a subset of 61 fields in 51 galaxies from 
the ACS Nearby Galaxy Survey Treasury (ANGST), a 
volume-limi ted sample of non-Lo cal Group galaxies out 
to ^ 4 Mpc ( Dalcanton et al. 2009 ) , and from the archival 
legacy program AINGRR R: Archival Nearby Galaxies: 
Reduce, Reuse, Recycle ( [Gilbert et al. 20111. Based 
on their HST stellar photometry, we have verified that 
these target fields have (1) sufficient numbers of stars to 
contain hundreds of candidate AGB stars within a sin- 



gle WFC3/IR FOV; (2) uncrowded stellar photometry 
of sufficient quality to provide useful constraints on the 
star formation history and metallicity distributiorp] (3) 
indications from the range of red giant branch colors that 
the stellar populations host a non-negligible fraction of 
stars whose metallicities fall outside of the range of well- 
studied LMC/SMC metallicities; and (4) a sufficiently 
broad range of star formation histories to allow sampling 
of evolving stars of different masses. Our target list in- 
cluded multiple pointings in the more massive systems, 
whenever the optical data suggested that there are signif- 
icant variations in stellar age and metallicity at different 
locations within the galaxy. 

During Cycle 17 we obtained observations for 26 of the 
61 possible SNAP targets. Two galaxies (Holmberg II 
and NGC 2403) had multiple pointings (2 and 3 point- 
ings, respectively), isolating regions with different star 
formation histories and metallicities. The properties of 
the 23 observed galaxies are listed in Table [T] Galaxy 
names, positions, apparent blue magnitudes (St), di- 
ameters, morphological T-types, and Hi line widths 
(W^so) have been ad opted from the primary name in the 
Karachentsev et al. (2004) Catalog of Neighboring Galax- 
ies! Distance moduir"are b ased on the F81AW TRGB 
from Dalcanton et al. I (12009) for most galaxies, and from 
Karachentsev^ al. (20^iorNGC_7793. Group mem- 
bershi ps are trom jKarachentsev ( 2005 ) or Tully et al 



(2006|). Foreground extinctions (Ay) are from ISchlegel 



et al.| ( 1998 1, as reported by the on-line Galactic Dust 
Extinction Service at the Infrared Science ArchivcP^ 

2.2. Observations 

Observations for this program were carried out during 
Cycle 17 as SNAP-11719. Some fraction of these were 
taken during the period of instrument commissioning af- 
ter the Hubble repair mission. All targets were observed 
in both FllOW^ and F160I^ with HST's WFC3 IR chan- 
nel. These filters offer the greatest depth in a given ex- 
posure time, and thus are likely to become the workhorse 
filters for the WFC3/IR camera. Observations were car- 
ried out in a 3-point "WFC3-IR-DITHER-LINE" pat- 
tern, with exposures in both i^llOH^ and F160H^ taken 
at each pointing. We adopted the STEP50 exposure se- 
quence, which accommodates data with a large dynamic 
range by using both short and long non-destructive reads. 
We used NSAMP = 9 in i^llOW^ and NSAMP = 11 in 
F160M^, giving total exposure times of 597.7 s and 897.7 s 
respectively. The exposure sequence was interleaved to 
allow buffers to dump during the exposures, with no la- 
tency. To maximize the schedulability of our observa- 
tions, we did not specify a roll-angle orientation con- 
straint. 

Properties of the actual observations can be found in 
Table ^ where we include both the name of the galaxy 
from Table [T] and the name of the specific target within 
the galaxy. Target names were chosen to match the 
names of existing optical data sets at the same position; 

Galaxies beyond 4 Mpc are sufficiently crowded that the 
deep optical photometry needed for deriving SFHs is not possible. 
Galaxies closer than 1 Mpc are sufficiently extended that few AGB 
stars actually fall in a typical FOV and contamination from Galac- 
tic foreground sources is severe, making this study more effective 
in more distant galaxies. 

http : / /irsa. ipac . caltech.edu/applicatlons/DUST/ 
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these associated data sets are also listed in Table [H The 
centers of the WFC3 field-of-view (FOVs) were chosen 
to maximize the overlap with the optical data. The lo- 
cations of the WFC3/IR footprints are shown in red in 
Figure [T] superimposed on Digitized Sky Survey images 
of each galaxy. Blue regions show the locations of the as- 
sociated optical imaging from ACS or WFPC2. Because 
we did not specify an orientation, occasionally a small 
corner of the WFC3/IR FOV fell off the area covered by 
the optical data. This mismatch occurred in cases where 
the optical FOV was not centered on the most desirable 
part of the galaxy, and the scheduled roll angle happened 
to be unfavorable; however, the fraction of the WFC3/IR 
FOV lacking optical coverage is less than 10% in almost 
all cases. 

In Figure [2] we show false color images of the obser- 
vations. There are a number of features to note. First, 
in many cases it is clear that our short exposures have 
already reached the "crowding limit" , where stars are suf- 
ficiently close on the sky that fainter stars could not be 
detected, even with longer exposures; this effect can also 
be seen as spatially-varying depth in the CMDs. Second, 
background galaxies are far more prevalent than in opti- 
cal HST images, suggesting that unresolved background 
galaxies are likely to be a more significant contaminant 
of the NIR CMDs. Third, the images are marked with 
a number of circles, located at positions of "IR blobs", 
that are flagged as bad data by the WFC3/IR pipeline 
(see the WFC3-2010-06 ISR by N. Pirzkal). Fourth, 
extended line emission is visible in i^llOiy images of 
young star forming regions, as can be seen in the close- 
up of a star forming region in IC2574-SGS, shown in 
Figure [3] The line emission in this filter is dominated by 
the [Siii]9069A,9532A doublet, with additional contribu- 
tions from Hel at 10830A and 10833A (from 2 photon 
emission) and Paschen /3 at 12818A. 

Finally, and perhaps most concerning, a number of the 
images show significant scattered light, particularly in 
the lower left region of the image (see UGC4305-1 or 
NGC2403-HALO-6, for example). We now discuss the 
origin of this scattered light component. 

2.2.1. Scattered Light 

Figure |4] shows the series of 6 combined *FLT images 
taken during our 1 orbit visit, after all detected sources 
have been fit and subtracted from the images. The se- 
ries shows images in both _F110VF (exposures 1, 3, & 
5 in the sequence), and Fl&QW (exposures 2, 4, & 6). 
The images show two types of structures in the "sky" 
background. The first is persistent structure due to un- 
resolved stars in the image; this component does not 
change throughout the visit, but varies from galaxy-to- 
galaxy. The second structure (falling typically on the 
left-hand side of the image, and peaking in the bottom 
left corner), varies rapidly in amplitude during a single 
orbit but has a consistent shape when present. The pat- 
tern of the time variation is not consistent, however, such 
that the brightness of the scattered light feature peaks at 
different times during an orbit. This variability leads to 
an apparent variation of the color of the scattered light; if 
the scattered light peaks during a FllOW^ exposure, the 
scattered light can appear blue, but if it peaks during an 
i^l60I^ exposure, the light can appear red. 



The rapid variation in the amplitude of the scattered 
light places strong constraints on its origin, as there are 
few telescope pointing attributes that vary significantly 
during an orbit. The angles between the principal VI 
axis and the Sun or the moon (FITS header keywords 
SUNANGLE and MOONANGL) vary little during an orbit, nor 
does the amount of zodiacal light. The only quantities 
which do vary on that timescale are the angle between 
the Sun and the Earth's limb (keyword SUN_ALT) and, 
more directly, the angle between the spacecraft pointing 
and the Earth's limb (keyword LQS_LIMB in the *.JIT 
or *.JIF files). We thank Ben Weiner (private commu- 
nication) for pointing out this possibility, based on his 
analysis of similar effects seen in WFC3/IR grism datgp^ 

To examine the variation of the scattered light as a 
function of the angle to the Earth's limb, we analyzed 
the sky levels in each star-subtracted * . FLT image. We 
assume that there are three contributors to the star- 
subtracted "sky" image: (1) a mean uniform background; 
(2) a spatially-structured, time-invariant background due 
to unresolved stars; and (3) a spatially-structured, time- 
variable background due to scattered light, which peaks 
in the lower left quadrant. To assess these different com- 
ponents, we recorded the mean level in the upper right 
quadrant of each image, which is assumed to be free of 
the scattered light feature, and in the lower left region 
dominated by the feature ([0:308,0:507], in image coordi- 
nates). Note that we do not expect these two sky levels 
to be identical, even in the absence of the scattered light 
feature, because of the contribution of light from unre- 
solved stars. However, they both should share a common 
"floor" , representative of the mean sky background dur- 
ing the observation. 

In the left hand panel of Figure [5] we show the mean 
sky level in the upper right quadrant of each exposure, 
plotted as a function of the limb angle between the space- 
craft and the earth, for the FllQW and FIQQW filters, 
respectively. If the earth's limb was dark at the time of 
observation, the sky level is plotted with an open circle. 
Observations taken within a single orbit are connected 
with a line. 

From the left hand side of Figure [5] it is clear that 
the sky brightness depends strongly on the limb angle 
for the FWQW filter, at least when the earth's limb is 
bright. Below angles of 40°, the mean sky level increases 
by roughly a factor of two or more. We expect that this 
correlation would be significantly tighter if this test were 
repeated for empty fields, where there would be little 
contribution from unresolved sources to the measurement 
of the sky. 

In the right hand panel of Figure [5] we attempt to cap- 
ture the amplitude of the scattered light feature visible in 
the first two panels of Figure|4j Unfortunately, we cannot 
simply track the variation ofthe sky level in that region, 
since we expect some limb angle dependent variation in 
the overall sky level, as seen in the left hand plot of Fig- 
ure [Sj Instead, we look at the limb angle dependence of 
the excess light in the lower left, compared to the upper 
right (Gradient = ^'\^y lower left - ^^Y upper right)- We as- 
sume that when the scattered light feature is absent, the 
difference in the sky level between the upper right and 

" See also ISR-ACS 2003-05 by Biretta et al for an analysis of 
ACS background light as a function of limb angle. 
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lower left regions reflects the structure in the sky due to 
unresolved sources. This difference should be constant 
with time, no matter what the mean sky background 
level is. For each set of observations (3 per filter in a sin- 
gle orbit), we then look at the "excess" in the difference 
between the lower left and the upper right, compared to 
the minimum during the orbit. We then scale the result- 
ing excess by the minimum sky level in the four quad- 
rants of the image, to give an indication of the fractional 
strength of the feature, compared to the typical sky level 
when unresolved sources are not present. 

The right hand side of Figure |5] shows that for most 
of the exposures, the relative sky brightness of the lower 
left and upper right varies by less than 1% between ex- 
posures, as long at the angle to the bright earth limb 
is greater than 30°. If the angle is less than 30°, how- 
ever, a strong scattered light feature is produced in the 
majority of cases. The onset is quite sharp, and is 
present in both filters, unlike the increase in the uni- 
form sky brightness with limb angle, which is gradual and 
much stronger in FllOW. Based on this analysis, the 
galaxies with the most prominent scattered light features 
(>4% in either fiher) are UGC4305-1, IGC2574-SGS, 
DD071, SN-NGC2403-PR, NGC2403-DEEP, NGG2403- 
HALO-6, and KDG73, in order of strongest to weakest. 
This ranking agrees with the visual impression seen in 
the color images in Figure [2j 

For our subsequent analysis, we make no additional 
attempt to compensate for the structure of the scat- 
tered light feature in these images. Since our photom- 
etry uses a very localized estimate of the sky, our results 
should be insensitive to this feature, except through hav- 
ing to accept more noise in the affected region. Given 
the brightness of the stellar sources, including observa- 
tions with low limb angle observations produces more 
benefits through increased integration time than are lost 
through increased noise levels. For many other possi- 
ble WFC3/IR programs, however, this may not be the 
case, and observers should consider specifying a mini- 
mum earth limb angle for their observations. 

There are two additional complications that may re- 
sult from the presence of rapidly varying scattered light, 
however. The first complication is that some reads af- 
fected by scattered light may have pixels that are erro- 
neously flagged as cosmic rays, if flux in the scattered 
light feature appears as an outlier compared to the large 
number of unaffected reads; however, we do not see in 
any obvious increase in the number of masked cosmic 
rays in the limb brightened regions in Figure |4j The sec- 
ond complication results from the fact that the final IR 
image is produced by fitting a linear function to a series 
of non-destructive reads taken during the course of the 
exposure. However, while the flux in a pixel from a star 
will rise linearly during the exposure, the flux from the 
background will not rise linearly when the angle to the 
limb or the limb's brightness changes during an obser- 
vation. The assumption of a linear fit may therefore be 
invalid when the flux in the background is comparable to 
the flux of the source in a given pixel. One should fur- 
ther note that this effect could be present even when the 
strong scattered light feature is not obvious, since, as one 
sees in Figure [5j the uniform sky background can vary 
strongly during the observation, particularly in FllOW . 
Unfortunately, it is not straightforward for us to evalu- 



ate this effect in our data. First, it is difficult to sep- 
arate the true time-variable sky background, from the 
spatially-variable time-constant background from unre- 
solved stars. Second, we cannot use spatial variations 
in the color magnitude diagram to identify systematic 
errors in the photometry associated with the scattered 
light feature, since galaxies typically exhibit significant 
gradients in the their stellar populations. 

2.3. Data Reduction & Photometry 
Photomet rv was perfor med using the DOLPHOT 



package[^ polphin] [2000| ) , including a new WFC3 
specific module generated in support of this HST pro 
gram. The core functionality of DOLPHOT remains as 
in previous versions (i.e., simultaneous PSF fitting across 
a stack of images aligned to sub-pixel accuracy), but 
the WFC3 update includes several new capabilities: a 
preprocessing routine that uses the data quality exten- 
sions to mask the image and then applies the appro- 
priate WFC3 pixel area maps, a WFC3 point spread 
function (PSF) library, WFC3 distortion corrections, and 
new zero points for WFC3. Pixel area maps, distortion 
terms, zero points, and encircled energy corrections were 
obtained from the WFC3 documentatioiC3 The PSF 
library was computed using the preliminary version of 
Tiny Tim 7.0 with WFC3 supporl[^ note that the cur- 
rent version of the WFC3/IR Tiny Tim PSF library is 
based on pre-flight data. For WFC3/IR, PSFs are com- 
puted with 10x10 subsampling and a 3" radius, for 64 
regions on the chip. DOLPHOT also adopts small aper- 
ture corrections calculated from uncrowded stars in the 
image; these corrections were always <0.01 mag. 

All PSF photometry was performed by minimizing 
residuals on all individual fit exposures, and then com- 
bined into final magnitudes for each star in each fil- 
ter. Non-detections are considered to be stars with 
signal-to-noise ratios of less than 4, for a given filter. 
Quality parameters are computed and included in the 
combined photometry as well. Two catalogs were then 
produced for each field using the DOLPHOT signal- 
to-noise, sharpness, and crowding parameters. The 
first "* . St" catalog contains all sources detected in at 
least one band with S/N>4, where the sharpness of the 
source does not exclude the possibility that it is stellar 
(sharpness^ <0.1). The second "*.gst" catalog only 
contains the highest quality photometry (including de- 
tections with S/N>4 in both filters, (sharpnessi^novK 
+ sharpnessi?i@ovi/)^ <0.12, and crowdingpnovy + 
crowdingf^igovF <0.48). The *.gst catalog typically 
contains ~'50±10% of the sources from the * . st catalog. 
These final parameter cuts were chosen to produce the 
cleanest CMD features (few stars outside of the known 
features) without culling large numbers of sources from 
the features themselves; the culled sources are almost 
entirely valid stellar detections, but have more uncertain 
photometry. The number of stars detected varies be- 
tween 6,000 and 40,000 per pointing, with a median of 



http : / /purcell . as . arizona . edu/dolphot/ 

Wb'CS PUb's from release of 15-iNov-2009; WFC3 zero points 
and eorrections to infinite aperture from J. Kalirai 2009 (ISR 
2009-30); WFC3 IDC tables use for distortion correction from 
t20100519_ir_idc.fits 

http:/ /www.stecf.org'/instruments/TinyTim/ 
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^-^10,000 (Table [2]). We also include the maximum and 
minimum stellar surface density of the * . st detections, 
calculated within 100 subregions defined by a 10x10 grid 
in the imager^ observations with a higher stellar surface 
density will be more affected by crowding, and observa- 
tions with larger differences between the maximum and 
minimum surface density will have stronger spatial vari- 
ations in the depth of the CMD and the accuracy of the 
photometry. 

Finally, we assessed the biases in our photometry using 
artificial star tests. False stars were placed into the image 
stack 100,000 times, such that the artificial stars sampled 
the full range of colors, magnitudes, and locations of our 
real photometry. We recovered the artificial stars using 
identical photometry routines and post-processing cuts 
as those applied to the data. 

Table [2] includes the resulting "50% completeness" 
magnitude for each filter, as determined by the limit 
where 50% of artificial stars inserted into the image 
are successfully recovered by our photometry (described 
below) in both filters (as in the * . gst catalog) . The 
50% completeness limits of our least crowded fields 
are -26.2 mag in FllOW and --25. 1 mag in F160W; 
these limits were significantly fainter than originally ex- 
pected for photon-limited sources, based on the pre-flight 
WFC3/IR Exposure Time Calculator. Note, however, 
that this depth is frequently position dependent, as can 
be seen from the position-dependent CMDs in the right 
hand panels of Figure [2j In crowded regions, the prox- 
imity of adjacent stars compromises the recovery of faint 
stars, producing brighter limiting magnitudes than one 
would expect based on photon counting statistics alone. 
The field for SN-NGC2403-PR has been particularly af- 
fected by crowding, and has a 50% completeness limit 
that is more than 2 magnitudes brighter than the me- 
dian of the sample. 

Figure [6] shows the magnitude uncertainties reported 
by DOLPHOT for the FllOW and F160T^ filters, for the 
full stellar "* . st" catalog. The small number of points 
that scatter to larger uncertainties typically come either 
from regions that were not covered in all three dithered 
exposures, or that have been fiagged as being poorly pho- 
tometered. The number of stars with higher than average 
uncertainties are reduced in the culled "*.gst" catalog, 
which produces somewhat sharper features in the CMD, 
at the expense of a modest drop in completeness. The 
effect of this culling can be seen in Figure [Tj where we 
plot the CMDs for the * . st and the * . gst photometry 
catalogs for one of our fields. We will describe the many 
features in these CMDs in Section |3] below. 

In Figure [8] we show the distribution of F160Vl^ mag- 
nitude errors for artificial stars inserted into our images, 
in bins of 1 magnitude, for a crowded and an uncrowded 
field (left and right plots, respectively). At the extremes, 
the distributions show a tail towards brighter recovered 
output magnitudes, due to the blending of faint unde- 
tected stars with brighter stars. In the median, how- 
ever, there is very little bias in the recovered magnitude 
(< 0.02 mag in the median for the faintest bin, which is 
nearly an order of magnitude smaller than the photo- 

Empirically, the data in Table [2] indicates that one cannot 
expect to detect more than ~5.5 stars per square arcsecond in an 
WFC3/IR frame. 



metric uncertainty at the same magnitude). Recovered 
colors are even less biased, as crowding typically offsets 
the magnitudes of both filters in similar directions. The 
amplitudes of these effects are of comparable amplitude 
in FllOW, but are not quite identical due to the different 
level of crowding and sky brightness in the two filters. 

2.4. Matching NIR and Optical Catalogs 

We match stellar catalogs from the distortion corrected 
optical and WFC3 images as follows. We calculate the 
astrometric transformation by first requiring —150 stars 
that are bright in both the optical and NIR data sets, and 
that spatially span the entire overlap region between the 
WFC3 and ACS (or WFPC2) images. We produce a list 
of candidate alignment stars by first culling the optical 
and NIR star catalogues from DOLPHOT to only stars 
that are in the spatial overlap region. We then select 
all reasonably bright, red stars in each dataset (opti- 
cal color > 0.7 mags and F814W < 26 mags; IR color 
FllOW - F160W > 0.5 mags and F160W < 24 mags), 
producing lists of —1000 stars in the optical and in the 
NIR. We sort the resulting two lists by luminosity to 
preferentially select luminous stars. Rather than select- 
ing the 150 brightest stars, which are frequently spatially 
clustered and do not span the full chip, we insist that as 
long as we have more than one star to choose from, we 
will select stars with an average separation of — 5", en- 
suring broader areal coverage of the alignment stars. 

We derive the transformation between the optical and 
IR star lists by first visually identifying a roughly linear 
shift between the two coordinate systems. We then itera- 
tively calculate a final tra nsformation using th e "method 
of triangles" described in Valdes et al. (19951, as imple- 
mented in the routine MATCH\'^'^\hy Michael Richmond. 
First we run a linear fit between the two coordinate sys- 
tems and then we use that solution as a starting point 
for a quadratic solution. We find that a cubic solution 
is generally unnecessary for the transformation between 
the distortion corrected WFC3 and ACS images, but we 
do use it for the transformation of the two data sets of 
WFPC2 observations. We apply the transformation to 
the entire NIR dataset bringing it into the optical coor- 
dinate system used for the ANGST and ANGRRR data 
releases. Note that we make no attempt to tie the images 
to the global astrometric frame, due to the lack of appro- 
priate astrometric standards in the small WFC3 fields of 
view. 

Once the WFC3/IR catalog has been astrometrically 
aligned to the optical catalog, we match individual stars 
in the two catalogs. We consider stars to be a match 
if the angular separation between the stars is less than 
0.07" (i.e., -0.5 WFC3/IR pixel). Typically 90% of the 
stars in the NIR catalogue are well matched to a star in 
the optical catalogue. Unmatched stars typically fall in 
the ACS chip gap or in the diffraction spikes of optically 
saturated bright stars. 

2.5. Characterizing the Age of the Stellar Population 

The morphology of the CMDs that result from our NIR 
and optical photometry reflect the age and metallicity 
of the underlying stellar population. To provide an in- 
tial constraint of these parameters, we have analyzed the 

http : //spiff . rlt . edu/match/match-0 . S/match.html 
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star formation history (SFH) of these galaxies using op- 
tical CMP fitting, sini ilar t o the procedu r es des cribed in 



Williams et al.| ([2011 1 and |Weisz et al.| ( |2011[) but re- 
stricted to the optical data that overlap the WFC3/IR 
FOV. 

Specifically, we measured the star formation rate and 
metallicity as a function of stellar age, by fitting th e op 



tical CMDs using the software package MATCH (Dol- 
l^in 2002). We adopted magnitude cuts set to the 507c 
completeness limits, and then fit the CMDs using linear 
combinations of the stellar evolution models of Girardll 



et al.| (2002, with updates in |Marigo et al.||2008| and Gi 



rardi et al.^2008J , populated with stars iollowmg an IMF 



with a slope of -2.3 and a binary fraction of 0.35 (with 
random mass sampling) , convolved with the photometric 
error and completeness statistics derived from artificial 
star tests. We first fit the data assuming a single fore- 
ground reddening a nd distance, adopting v alues used in 
the ANGST survey (palcanton et al.|2009[ ) based on the 



Schlegel et al. (1998) Galactic dust maps and the mag- 
nitude of the i KGB, respectively. The best fit provides 
the relative contribution of stars of each age and metal- 
licity in each field, which is then converted into cumu- 
lative stellar mas s produced a s a fun ction of time. For 
consistency with Weisz et al. (2011), we restricted the 
metallicity evolution to be constant or increasing with 
time, for galaxies in common with the Weisz et al. (2011 1 



dwarf sample. This choice provides more accurate recent 
star formation histories, at the expense of introducing 
occasional biases against SF in the oldest bin for the few 
galaxies with very deep data 

To assess the uncertainty of the best fit, we ran ex- 
tensive Monte Carlo tests. These tests assess two types 
of uncertainties: (1) random errors due to Poisson sam- 
pling of the CMD and errors in photometry; and (2) sys- 
tematic errors due to deficiencies in the stellar evolution 
models, or due to offsets in distance, reddening, and/or 
magnitude zero-points. The Poisson errors are accounted 
for by generating artificial CMDs from the best-fit con- 
volved model 100 times, and refitting the resulting CMD. 
While fitting each of these realizations, the systematic er- 
rors are assessed by introducing small random shifts in 
log(roff) and Mboi- The random values were drawn from 
a Gaussian distribution with a width of 0.03 in log(Tcff) 
and 0.41 in Afboi- The size of these shifts are set by dif- 
ferences between models in the literature, and therefore 
serve as a proxy of the effects of our particular choice 
of stellar evolution models. Our final uncertainties are 
the 68% confidence intervals of the results of all of our 
Monte Carlo test fits. These total uncertainties are used 
as the error bars in all subsequent plots and analysis. 
Note however that adjacent time bins in the SFH are 
covariant to some degree, such that when one time bin 
fluctuates high, the adjacent bins fluctuate low. As a 
result, the plots of the cumulative SFH offer a truer rep- 
resentation of the uncertainties. 

The resulting SFH and cumulative age distributions 
are shown in the bottom panels of Figure |9] The lower 
left panel shows the SFR as a function of time, and the 
lower right hand panel shows the cumulative star forma- 
tion history, measured from the present back to 14 Gyr 
ago. A detailed discussion of the SFH analysis of the 
NIR CMDs will be discussed in an upcoming paper. 



3. COLOR MAGNITUDE DIAGRAMS 

In this section we describe the properties of stellar pop- 
ulations in the FllOW + F160W filter set using both 
models (Section 3.1) and observations (Section |3.2|. We 



discuss the qualitative dependence of these properties on 
the age of the stellar population, as derived from optical 
CMDs in Section [231 



3.1. Overview of Key Features in Model CMDs 

To elucidate interpretation of the CMDs, in Figure [TO] 
we show simulated NIR color magnitude diagrams, color- 
coded either by age or stellar mass (for a constant star 
formation rate at a fixed metallicity of [Fe/H]=-1.45), or 
by metallicity (for a burst of star formation between 8.9 
and 11.2 Gyr). These simulations include the updated 
models for AGB ma ss-loss at low metallicity described in 
Girardi et al. (20101, and use the artifical stars from the 
M81-DFFP field to calculate typical photometric uncer- 
tainties and bias es. Comparable figures for optical CMDs 



can be found in Dalcanton et al 
In both the optical and the IN 



mm. 

K, the most prominent 



feature is the red giant branch (RGB), found at colors of 
FllOW-FieOW^O.S in the NIR and F606VF-F8141^- 
1 or FA75W - F8UW - 2 in the optical. This color 
depends on metallicity (Sect ion [6|), and beco mes bluer 
when the metallicity is low (Aaronson et al. 1978) . At 
low metallicities, the optical KGB also becomes more 
vertical and exhibits less curvature. In the NIR, however, 
the slope of the RGB is ne arly vertical, with on ly modest 



variations with metallicity (Ferraro et al. 2000[ ). The ages 



of stars in the RGB can span a wide range (> 1 Gyr) 
For the galaxies in o ur sample, early star formation is 
partic ularly vigorous (Weisz et al. 2011 Williams et al.| 
20111, which will weight the population of RGB stars 



towards older ages, compared to the constant SFR shown 
in Figure [TOj 

One of the most prominent features along the RGB 
is the "red clump" , typically found at 3-4 magnitudes 
fainter than the TRGB in the optical. The stars in the 
red clump are burning Helium in their cores in the same 
way as horizontal branch stars, but appear red due to 
either their young age s (large envelope mass) or their 
high metallicities (see Girardi & Salarisj (2001 ); Salaris 
f2002); CastcUani ct al. (2000) for theo retical models 
and Ivanov & Borissova, (.2002J ; ^Grocholski fc Sarajedini 



( 2002| ; [Valenti et ar i 2od4|) tor observational constraints 
from globular clustersJT Unfortunately, our NIR data 
does not reach this information-rich feature, due to the 
shortness of our NIR exposures and the bright crowd- 
ing limit resulting from the larger WFC3/IR pixel scale. 
The absence of this feature reduces the utility of using 
our NIR CMDs to constrain the relative amounts of an- 
cient and intermediate age SF. 

Ancient star formation also produces a population of 
low mass thermally pulsing AGB stars (TP-AGB; < 
1.5 M0), found just above the TRG^at Mf814W~-4. 
At younger ages, TP-AGB stars are more massive and 
have much brighter magnitude s (iMarigo et al.|2008 ) . As 
we discuss below in Section |3.2.1[ the colors of AGB 

AGB stars are also present in the same region occupied by the 
RGB, but are greatly outnumbered. They cannot be distinguished 
cleanly as a separate population fainter than the tip of the red 
giant branch at MFgi4w~— 4. 
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stars vary little in the NIR, and thus produce nearly 
vertical sequences for all bu t the smaller population of 
red extreme AG B stars (e.g. 'Ni kolaev fc W einberg"2000| 
GuUieuszik et al. 2008; Melbour ne et al.||2 010a; Davidge 



2om 

Figure [To] shows that younger stars are expected to 
dominate the stellar populations blueward and bright- 
ward of the RGB. The youngest stars are main sequence 
stars, found at the bluest edge of the CMD. While the 
main sequence is quite well defined in the optical, it is not 
nearly as distinct in the NIR. This difference must result 
in part from the fact that only the most luminous high 
mass main sequence stars are detectable in these NIR 
observations. In the optical CMDs, we typically detect 
main sequence stars with masses of 4.5 M0 and above, 
whereas in the NIR, we can typically only detect main se- 
quence stars with > 20 Mq (assuming a magnitude limit 
of 25 mag in F160W versus 28 mag in F814:W, for a me- 
dian distance of 3.5 Mpc and metallicity of O.IZq). Such 
stars are typically rare, and thus the detectable main se- 
quence will not be sufficiently well-populated to appear 
as a distinct sequence unless there has been ample very 
recent star formation. 

Complicating the clear detection of the main sequence 
in some cases is the presence of slightly older blue he- 
lium core burning (BHeB) stars. These are evolving 
post-main sequence 5 — 2OM0 stars, with ages of typi- 
cally 10 — 500 Myr. These stars, along with their red 
counterparts (RHeB stars), are formed following the ex- 
haustion of core hydrogen burning, when the core rapidly 
collapses and the stars' effective temperature decrease^^ 
The core soon reaches an equilibrium, when core helium 
burning and inner shell hydrogen burning commence, ex- 
panding the redder outer layers. A star at this phase is 
known as a red helium burning star. As helium in the 
core is converted into carbon, the RHeBs become visibly 
hotter and enter the blue helium burn ing (BHeB) phase 
of evolution (e.g., Langer fc Maeder[T9 95 ) . Although the 
BHeB phase occupies a lar ge traction of the total core he- 
lium burning lifetime (e.g.,lBertelli et al.|1994), stars are 
capable of alternating between HHcB and RHeB stages, 
with the precise lifetimes determined by the intricate re- 
lationships betwee n underlying physical parameters (e.g., 
Chiosi et al.|[l992| . 

'I'he red and blue core helium burning sequences are 
most obvious at optical wavelengths. They appear as 
two sequences emerging brightward of the red clump, 
where low mass < 2 Mq core helium burning stars are 
found. The BHeB sequence emerges diagonally, heads 
to bluer colors at higher luminosities, then becomes ver- 
tical and parallel to the main sequence at high lumi- 
nosities and high stellar masses The RHeB sequence 
emerges nearly vertically from the red clump, and ex- 
tends steadily brightwards, at somewhat bluer optical 
colors than the RGB. Along these sequences, stars of 
a given mass appear at a single luminosity, so that the 
number of stars along the sequences indicate the numbers 
of evolving stars of different masses. This connection al- 
lows the star formation rate between ~5 and ~400 Myr 



There is some overlap in the literature between the brightest 
stars in the HcB sequences and classical red and blue supergiants. 
However, the former extends to much lower masses and luminosi- 
ties. 



to be essentially read off the se quence, as was pioneered 
by Dohrn-Palmer et al. ( 2002 1, and more recently dis- 
cussed in McCjuinn et al.| ( 2011 ) . More massive, and thus 
younger, stars appear at very bright magnitudes along 
both the BHeB and RHeB sequences. However, we note 
that the models for these stars are currently uncertain, 
and are known not to produce the correct colors or rela- 
tive number of stars on the red and blue sequences at low 
metallicities (Dohm-Palmer & Skillman 2002 McQuinn 
eraL]|20TI| . 



3.2. Observed Properties of NIR CMDs 

We now turn to the observed NIR and optical CMDs 
for the fields covered by our WFC3/IR observations, 
shown in the top panels of Figure [9] In the upper left 
panel of Figure |9] we show the NIR (JMD of the cleaned 
* . gst photometry catalogs, uncorrected for foreground 
extinction or reddening. In the adjacent panel, we show 
the highest quality optical CMD that is available from 
overlapping archival imaging. The stars in the optical 
CMD have been restricted to those that overlap the area 
covered by the WFC3/IR FOV. These observed CMDs 
can be compared with the models in Figure [T0| to under- 
stand which phases of stellar evolution are populating 
various features in the CMD. We now discuss the princi- 
pal features of these CMDs. 

3.2.1. Old and Intermediate Age Populations: Ttie AGB & 

RGB 

In a broad qualitative sense, the NIR CMDs in Fig- 
ure [9] show the features expected from the stellar mod- 
els, particularly at old ages. In fields dominated by old 
stellar populations (as judged from the optical CMD; see 
KDG 63, or NGC 404 for examples), the NIR CMDs 
are quite simple. They show roughly 3 magnitudes of 
a well-defined upper red giant branch terminating at 
-^^Fieow ~ —5.7. The red giant branch is more vertical 
and exhibits less curvature than in the optical (see M81- 
DEEP for an example), as expected from Figure 10 The 
color at the tip of the RGB varies from between CTy and 
1.1, as we discuss in more detail in Section|4]below. This 
range of colors agrees well with the range expected from 
the models in Figure [TO] Unfortunately, the data are not 
deep enough to reveal the red clump in the NIR; this 
feature is expected to appear ~4 mag below the TRGB. 

Galaxies dominated by old stars also host a sparse pop- 
ulation of TP-AGB stars located above the tip of the 
RGB (AfFi6ow<-5.7). The NIR population of TP-AGB 
stars falls in a vertical sequence spanning a narrow range 
of color, in contrast to the broad "fan" of TP-AGB stars 
seen in the optical. At optical wavelengths, the TP-AGB 
stars span a wide range of colors, due to a combination 
of dust in their circumstellar envelopes, large variations 
in molecular line spectra with photospheric temperature, 
the distance between the op tical and the NIR peak of an 
AGB's star bolometric flux (jFrogel et al. 1990| , and long- 
period variability (which aflects all passbands, but only 
produces significant color spreads in the optical). 

The NIR TP-AGB sequence typically extends < 1 mag- 
nitude in F160W above the TRGB. However, the TP- 
AGB reaches even brighter magnitudes in populations 
with more recent intermediate age SF (see NGC 3741, 
for an example). The existence of bright, massive TP- 
AGB stars at younger stellar ages is predicted by the 
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m odels as well (see F igure 10 and luminosity functions 
in |01sen etal] ( |2006| ). 

At lower luminosities in the TP-AGB sequence, CMDs 
sometimes show a noticeable "gap" between the bright- 
est RGB stars and the faintest TP-AGB stars above the 
TRGB (see NGC2403-HALO-6). This localized drop in 
the luminosity function seems to hav e a counterpart in 
evolu tionary tracks of TP-AGB stars ( Marigo & Girardi 
20071. The minimum likely reflects tne steeper rate of 
brightening (— dMboi/di) that characterises the initial 
stages of the TP-AGB phase, when the first thermal 
pulses have still not reached the full-amplitude regime 
and TP-AGB stars are expected to be fainter than pre- 
dicted by the core mass-luminosity relation. As a result, 
the brighter late stages of TP-AGB evolution are slower, 
such that luminosity bins appear relatively more pop- 
ulated than the fainter early st ages. Simulated CMDs 
based on the [Marigo fc Girardij ( |2007) TP-AGB models 
do show a local minimum in the luminosity function right 
above the TRGB, as expected. However, the ga p is not 
apparent in the more recent models desc rib ed in [Girardi 
et al. (20101) (lower right panel of Figure 10). The obser- 
vation ot this features opens the possibility that it could 
help to calibrate the models at the early stages of the 
TP-AGB. 

The mean color at the base of the TP-AGB sequence is 
comparable to the mean color at the tip of the RGB, such 
that the TP-AGB sequence appears to emerge vertically 
from the RGB. However, the cumulative distributions 
of colors immediately above and below the TRGI^^ sug- 
gests that the AGB is slightly redder by 0.02-0.04 magni- 
tudes for galaxies where there is minimal contamination 
from RHcB stars (e.g., KKH37, HS117, DD082, DD078, 
N2976). Visually, the TP-AGB sequence appears to span 
a smaller range of color than the RGB in some cases. 
However, there is no statistical evidence that this is the 
case, as both the RGB and TP-AGB have statistically 
indistiguishible widths, when the fewer number of stars 
in the AGB is taken into account. 

We note that our sample does not immediately show 
a large obvious population of carbon-rich AGB stars. In 
the corresponding J + H ground-based NIR filter set, 
such stars dominate a roughly diagonal sequence that 
starts near Mh ~ — 7 and extends to redder colors and 
fainter magnitudes. To illustrate this sequence, we show 
empirical data in Figure 11 using J — H vs H CMD for 
2MASS observations of tEe the LMC. This fiher com- 
bination is the closest match to the FllOW + F160W 
filter set used in our observation^^ In the LMC data 
the sequence of carbon-rich and/or extreme AGB stars 
is roughly horizontal at Mh ~ —7.3, and extends red- 
ward from the sequence of more numerous oxygen-rich 
AGB stars. There are a few cases where a compara- 
ble red sequence of extreme AGB stars is hinted at in 
our data (e.g., UGC4305, NGC4163, NGC2403-HALO-6, 

21 We define the RGB stars "below" the TRGB as being the 
0.25 magnitude interval starting 0.05 mag fainter than the TRGB, 
and the AGB stars "above" the TRGB as being the 0.5 magnitude 
interval starting 0.05 mag brighter than the TRGB. This selection 
produces roughly equal numbers of stars in each subsample. 

22 -pov reference, at J — _ff = 1, magnitudes in F160W are 0.28 
magnitudes fainter than in H, and colors in _F110Vl^ — F160W 
are 0.08 magnitudes redder than in J — H, as shown below in 
equations [2] & |3] 



NGC0300-WIDE1, DD082). However, in none of these 
is the sequence as noticeable as in the LMC. 

The absence of a strong carbon star tail is initially 
suprising given that our sample of galaxies is dominated 
by lower metallicity galaxies than the LMC. The dredge- 
up of carbon is believed to be more efficient at lower 
metallicities, both because of more favourable conditions 



for the occurence of strong helium shell flashes ( Karakas 
et al. [2002[ Stanchffe 2006), and because of tne lesser 
carbon dredge-up needed to reach a carbon-rich condi- 
tion (with C/0>1 at the surface). Thus, one expects 
carbon-rich stars to be more numerous in low metallicity 



eral studies (e. 
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2005 


Groenewe- 


gen|1999 


2007 


Boyer et al.|,2011). 



are significantly hotter than at solar metallicities, which 
inhibits the formation of the molecules that cause the 



spectral fe atures typical of carbon stars (Marigo & Gi 
rardi||200"7 ). The resulting red uced opacity in low met"al 



licity carbon-rich stars (e.g., [Marigo fc Aringer 20091 
may move the carbon stars back to bluer iNlR colors than 
their hig h metallicity counterparts (see, for instance. Fig- 
ure 7 of Marigo & Girardi 2007). Thus, while carbon 
stars may be more numerous in low metallicity systems, 
they may be less obvious outliers in NIR CMDs. The 
difficulty of identifying carbo n stars in the NIR has been 
demonstrated empirically by [Battinelli & Demers (20091 
for galaxies within the Local Group. 

We believe that the lack of obvious carbon stars is fur- 
ther exacerbted by use of the F 110W+F160W filter set. 



Simulations with TRILEGAL (Girardi et al. 



2005 G 



rardi & Marigo 20071 shown in I'igure T2 indicate that 
the majority of carbon-rich AGB stars occupy similar 
locations as oxygen-rich AGB stars in the WFC3/IR 
C MDs. These simul ations were p erformed using both 
the [Loidl et aLl ( [2001[ ) (top row) and [Aringer et al.| ( [2009[ ) 
(bottom row) synthetic spectra for C stars. Whereas car- 
bon stars are distinctly redder in J — H than the oxygen- 
rich AGB (left column), in FllOW - F160W they are (if 
anything) slightly bluer (middle column), particularly for 
the Ari nger et al.[ ( 2009j ) atmosphere models. The prop- 
erties of the carbon-rich AGB population in our sample 
will be explored in more detail in an upcoming paper. 

3.2.2. Young populations: The Importance of Red Core 
Helium Burning Stars 

In contrast to older stellar populations, young stel- 
lar populations have surprisingly complex NIR CMDs. 
Fields with the most active star formation show far 
more features than just the RGB and AGB (see targets 
IC2574-SGS (in 12574), NGC0300-WIDE1 (in N300), 
UGC4305-1 (in HoII), and UGC5139 (in Hoi)). We now 
identify and briefly discuss these various features. 

Star forming NIR CMDs host a vertical sequence at 
FllOW - F160W ~ on the blue side of the CMD, 
made up of both main sequence and BHeB stars. They 
also show a broad cloud of stars with colors interme- 
diate between the main sequence and the RGB. These 
stars are most likely to be BHeB stars and other evolv- 
ing massive stars. Note that only the most massive blue 
stars are detectable in the NIR. Such stars are rare, mak- 
ing it less likely that a "sequence" will be sufficiently 
well-populated to appear as a distinct feature (see Fig- 
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ure |10[ ), which reduces the clarity of the main sequence. 
In neids with on-going but low intensity star formation 
(e.g., targets NGC3741, NGC7793-HALO-6, NGC3077- 
PHOENIX), the "main sequence" manifests more as a 
blue edge to the cloud of points, rather than the narrow 
sequence seen in the optical. 

To elucidate the correspondence b etw een features in 
the optical and NIR CMDs, in Figure [13] we show CMDs 
for stars that have been matched between the NIR and 



but the slightly older clusters have large concentrations 
of extremely luminous RHeB stars. 

Thus, although much recent work has focused on the 
importance of the AGB phase to setting the NIR colors 



optical photometry catalogs ([2.4) for IC2574-SGS, after 
correcting for foreground extinction and distance. The 
stars have been color-coded by their likely evolutionary 
phase, as deduced from the optical. For this exercise, 
we use optical catalogs involving the F555 filter, which 
provides good separation of CMD features, particularly 
between main sequence and BHeB stars, and between 
RHeB and RGB stars. Stars that did not have high qual- 
ity positional matches (coincident within 0.07"), or that 
were faint in F81AW, are plotted in black. We note 
that the matching is unlikely to be perfect, and that 
some modest fraction of stars in the NIR catalog may 
have been incorrectly matched with stars in the opti- 
cal catalog. This mismatching would be most likely in 
the crowded star-forming regions in this particular field, 
where slight errors in the alignment could still leave stars 
within t he e rror radius for matching. 

Figure [T3| reveals a number of features. First, as ex- 
pected, the main sequence stars (in blue) crowd along 
the blue edge of the NIR CMD. However, the blue core 
helium burning stars (in green) also contribute the blue 
vertical sequence, as well as filling the region between the 
NIR vertical blue sequence and the RGB. Thus, the lack 
of a clear BHeB sequence in the NIR is due primarily 
to its merging with the MS. On the red side of the NIR 
CMD, we see the R GB (red) and AGB (cyan) features 
discussed in |3.2.1 along with a modest horizontal tail 
of likely carbon-rich AGB stars at Mfi6ow~— 7. In ad- 
dition to the AGB stars seen in the CMDs of older stellar 
populations there is a brighter plume of AGB stars ex- 
tending brightward of Afpieow ~ — 7.5 and AfF8i4W ^ — 5. 
This plume is likely due to the presence of more mas- 
sive stars on the AGB, as a result of the strong SF at 
t < 2 Gyr in the IC2574-SGS field. Note also that the 
AGB stars fall on a well-defined sequence in the NIR 
CMD, but spread over a wide swath of color in the opti- 



cal, as discussed above in Section 3.2.1 



The most remarkable feature in Figure [T3| is the strong 
sequence of luminous RHeB stars. The color of the se- 
quence is only slightly bluer than the AGB and RGB 
(by ^-^0.2 mag), and extends to far brighter magnitudes. 
Indeed, the red core Helium burning evolutionary phase 
is responsible for all of the most NIR luminous stars 
in IC2574. As can be seen in Figure |3) these lumi- 
nous RHeB stars are tightly clustered wTthin an active 
star forming region that features several young clusters 
spanning a ^ 100 Myr range of ages. Looking at the 
optical image (left), the youngest clusters in the region 
(< 10 Myr) are still embedded in luminous HII regions, 
while the slightly older clusters are somewhat fainter, 
more diffuse, and lack extended line emission, although 
they are still quite blue, due to large concentrations of 
luminous BHeB stars. However, in the NIR (right), the 
relative luminosities of these clusters are reversed, such 
that the youngest clusters have negligible NIR emission. 



and lum inositi es of young ste llar populations (Maraston 
et al.|2006; Henriques et al.|2'010j ) , it appears that RHeB 



stars can potentiall y b e equally important. We demon- 
strate this in Figure [l4j where we show the integrated lu- 
minosity of stars in different evolutionary phases. At the 
limit where our data becomes incomplete, RHeB stars are 
the dominant contribution, followed by AGB stars. The 
RGB, which is canonically assumed to dominate the NIR 
light, is only the third most important contributor to the 
luminosity at the completeness limit of our data. Note, 
however, that the fractional contribution of RGB stars 
will increase when the full range of stellar luminosities is 
considered; corrections fo r the missing stars are i ncluded 
in a companion paper by Melbourne et al. (2011 1, where 
the fraction of integrated light contributed by RHeB and 
AGB stars is quantified, showing that the RHeB can con- 
tribute up to 25% of the total light in i^l60T4^, and that 
current models underpredict the luminosity contribution 
of RHeB stars by up to a factor of 4. 

The contribution of RHeB stars to the total luminos- 
ity will be strongest when the SFR has been elevated 
between 25 and ~ 100 Myrs ago, which produces RHeB 
stars brighter than AfF8i4W ~ ~5. This timescale is the 
same as the one over which UV emission is expected to 
be significant, and thus surveys that select for galaxies 
with high UV flux may also be selecting for galaxies with 
a significant RHeB population. Failure to account for the 
luminous, somewhat blue RHeB population would lead 
one to overestimate the mass-to-light ratio in the NIR, 
and to infer lower metallicities from the broad-band col- 
ors. Unfortunately, the models of RHeB stars are even 
more uncertain than for AGB stars, and are known to 
produce erroneous colors and re lative num bers of BHeB 
and RHeB at some metallicities ( Langer fc Macdcr 19951 
Dohm-Palmer fc SkillmaHl|2002l [Gallart et al.||2005i |M^ 



Quinn et al.||2011j) 

Uur conclusions about the importance of RHeB stars 
are not significantly affected by uncertainties from fore- 
ground contamination. Unlike studies in the Magellanic 
Clouds, our fields cover small areas on the sky, which 
minimizes contributions from bright Milky Way stars. 
The lack of luminous contaminants can be seen empiri- 
cally in Figure [2] (presented below) , where we show the 
CMD of different subregions in the chip. We see essen- 
tially no luminous stars in the outer regions of small 
galaxies, where we are presumedly dominated by fore- 
ground stars and background galaxies (see HS117, for ex- 
ample). Pre dictions from the de fault TRILEGAL Milky 
Way model ( Girardi et al.|2005 ) also suggest that we ex- 
pect fewer than one contaminating star in each of the 
NIR CMDs plotted in Figure (91 

We also note that individual RHeB stars could eas- 
ily be confused as individual stellar clusters, in images 
of galaxies that are nearby, but that are not sufficiently 
close to resolve individual stars. Due to their high lumi- 
nosity and red colors, it may be difficult to distinguish 
single IR-luminous RHeB stars from older stellar clusters 
with larger numbers of fainter RHeB stars. 
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3.3. Color-Color Diagrams 

We can further explore the separation of different stel- 
lar evolutionary phases using the optical-NIR color-color 
diagram. We adopt the target IC1574-SGS as a test case, 
and use the matched catalogs from Figure [13] to gener- 
ate an optical-NIR color-color diagram, color-coded as in 
Figure 13 The resulting color-color diagram is shown in 



15 



Figure 

Figure 15] shows that the optical and NIR colors 
are highly correlated for colors bluer than F555W — 
i^814M^< 1.7. In this regime, the optical and NIR col- 
ors track each other extremely well, and are thus of lit- 
tle utility for separating phases of stellar evolution in 
data of this quality. At redder colors, however, there is 
noticably decreased correlation between the optical and 
NIR colors, driven almost entirely by the behavior of the 
AGB population (cyan points in Figure 15). Although 



the bluest AGB stars follow a narrow sequence of optical- 
NIR colors, the optically redder AGB stars are far less 
well-behaved. To first order, the NIR color saturates 
at i^llOM^ — F160W ~ 1 for a wide range of optical 
colors, with the e xception of a likely population of ex- 
treme AGB stars (Blum et al. 2006). Qualitatively, it 
appears that the narrower blue AgB sequence bends over 
to fixed NIR color with increasing optical color. Beyond 
i^555M^ — F814M^ > 1.7, the narrow sequence appears to 
become embedded in a larger swath of red AGB stars, 
for which there is little correlation between optical and 
NIR color, and wide dispersion. There is no qualitative 
evidence for the narrower sequence continuing redwards 
of - F8UW>2.2. 

We can further explore these different regimes of AGB 
behavior by qualitatively separating the AGB sequence 
into two regions and examining the behavior of the stars 
in the NIR CMD. Figure [IB] shows the adopted division 
of the AGB population in color-color space (upper left), 
and the resulting NIR CMD (upper right). Although 
the division into two classes is somewhat arbitrary, and 
not motivated by a particular choice of stellar model, 
the adopted separation does appear to break the AGB 
population into two luminosity classes. As shown in the 
histograms in the lower left panels, the redder, high- 
dispersion sub-population is a constant fraction of the 
AGB stars brighter than Mpieow < —6.8, it also lacks 
lower luminosity AGB stars, and makes up fewer than 
10% of the stars fainter than Mpigow > —6.2. In light of 
our selection criteria, this luminosity difference suggests 
that the higher luminosity, more massive AGB stars have 
systematically redder optical colors (lower right panel) 
and/or bluer NIR colors (with the exception of the ex- 
treme AGB stars redward of i^llOW^ - F160W > 1). 
Brighter than Mpiggw < —6.8, however, we see no noti- 
cible difference in the NIR luminosities or median NIR 
colors of the two sub-populations (beyond the larger dis- 
persion in the NIR color expected by the adopted division 
of the two populat ions in color-color space). 

Although Figure [16] suggests that there are systematic 
luminosity-dependent differences in the optical-NIR col- 
ors of AGB stars, we currently lack a physical model 
that would allow us to make better motivated choices 
for dividing the AGB population on the basis of optical- 
NIR colors. It is tempting to think that the division 
we have made is helping to isolate carbon stars, based 



on the behavior seen in Figure |12[ On the other hand, 
that same figure shows the large current uncertainties 
in the modeling at the relevant wavelengths. A more 
definitive approach to isolating and modeling carbon 
stars will likely require observations at somewhat longer 
wavelengths than are currently possible with WFC3's IR 
channel, or the addition of narrow-band imaging. 

3.4. Spatial variations 

In addition to the field-to-field variations of the NIR 
CMDs discussed above, we also find strong spatial varia- 
tions within individual fields. Such variation is naturally 
expected if young and intermediate age populations con- 
tribute significantly to the NIR CMD. In the right hand 
panels of Figure|2]we show the NIR CMDs in 16 subfields 
within each WFC3/IR frame, next to a color image of the 
same frame. 

These grids show a number of features. First, one can 
see the effects of crowding. In the most well-populated 
fields, the depth of the CMD is strongly dependent on 
position, such that the densest subregions have the shal- 
lowest depth. Second, one can see significant spatial vari- 
ations in the underlying stellar populations, particularly 
in fields that have significant amounts of recent star for- 
mation. Features due to the main sequence, RHeB, and 
AGB vary in strength and luminosity within individual 
fields. As an example, in NGC0300-WIDE1 the steUar 
populations from the left side of the frame host a promi- 
nent main sequence and luminous RHeB stars, whereas 
these populations are nearly absent from the lower right 
of the frame. Finally, CMD features due to young stel- 
lar populations typically appear with more clarity in the 
subregions than in the field as a whole (for example, note 
the sharp localization of the RHeB and BHeB in the up- 
per left quadrants of IC2754-SGS, or in the relative pop- 
ulation of the AGB and RGB across the field). 

4. TIP OF THE RED GIANT BRANCH IN THE NIR 

In the optical, the magnitude of the tip of the red gi- 
ant branch has become a widely used distanc e indicator 
for galaxie s with r esolved stellar populations ( Lee et al. 
1993. iSakai et al.|[T996. ,Mendez et al.,,2002l iK aractieni- 
sev et al.|2006||Makarov et al.|2006D. Its utility as a dis- 
ic[i( 



tance indicator rests on the relative insensitivity of the 
TRGB magnitude to age (for age> 3 Gyr) or metallicity 
(for [Fc/H]<— 0.5) in the / band, wher e the bolom etric 
magnitude of RGB stars typically peak ( Lee et aL|T993j) . 

In the NIR, however, the behavior of the i'RGB magm- 
tude is quite different. Theoretical models indicate that 
the magnitude of the TRGB should vary significantly 
with the properties of the un derlying stellar population 
(Figure 2 of Salaris & G irardi|200 5 ) , complicating the use 
of the NIR TRGB as a distaiice indicator. The variation 
in the NIR TRGB magnitude is known to be metallicity- 
dependent for sub-solar metallicities (< 0.5^0 ), as has 
been demonstrated conclusive ly in data for globula r clus- 



ters in the 2MASS filter set ( [Valenti et al.]2004| . One 
can therefore potentially correct tor the ^ 1 magnitude 
variation in the NIR TRGB absolute magnitude if the 
metallicity is known, allowing the TRGB to be used as a 
distance indicator in the NIR. 

Unfortunately, this basic approach may not be effec- 
tive in extragalactic systems, where the mean metallic- 
ity is uncertain, and the underlying stellar population 
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is more compl ex. As shown in Figure 10 of |Melbourne| 
et al. (2010b), the presence of intermediate age stars 
(1-2.5 (iyr) shifts the NIR TRGB to significantly fainter 
magnitudes. However, one could potentially diagnose the 
presence of this population through other means (for ex- 
ample, through a larger population of bright AGB stars) , 
and correct for it. The metallicity dependence could like- 
wise be corrected for as well; since metallicity directly 
affects the color of the RGB, it should be possible to em- 
pirically calibrate a relationship between the magnitude 
and color of the TRGB. 

In this section, we derive the magnitude and color of 
the TRGB for our sample. We assume that the distances 
in Table [I] (based on _F814W^ TRGB measurements) are 
correct, and then derive the absolute magnitude of the 
NIR TRGB. In what follows, all magnitudes and colors 
have been corrected for the foreground extinction given 
m Table |lj assuming t hat ApnnwMv = 0.33669 and 



^Fi60w/-5v = 0.20443 (jGirardi et al.[2008 , updated to 



included the latest in-hight calibrations tor the WFC3 
filter sets). 

As a first step, in Figure [IT] we plot the observed lu- 
minosity functions of the red extinction-corrected stars 
(dark black histogram), and their median color in mag- 
nitude bins (red line, plotted only for bins with more 
than 12 stars). Red stars are selected using a magnitude- 
dependent color cut to suppress the contribution of bluer 
young BHeB and MS stars; brighter than the TRGB, 
we keep all stars redward of FllOM^ — F160W^ = 0.5, 
and fainter than the TRGB, we keep all stars redward 
of the diagonal line connecting m = totrgb, FllOW — 
F160W = 0.5 with m = totrgb+4, F110W-F160W = 
—0.2. This cut has no appreciable effect on any feature 
visible in the luminosity function. The luminosity func- 
tions have not been corrected for incompleteness, and we 
therefore expect them to roll over at faint magnitudes due 
solely to observational e ffects; we carry ou t a ful l correc- 
tion for these effects in Melbourne et al. (2011), where 
we analyze the NIR luniinosity contributed by different 
evolutio nary phases. 

Figure [17] shows that there is a clear variation in the 
absolute magnitude of the TRGB, relative to a fiducial 
TRGB absolute magnitude of Mpieow — —5.7, plotted 
as a vertical dotted line. Visual inspection shows that 
redder RGBs typically terminate at brighter absolute 
magnitudes; for these presumably metal rich RGB stars, 
the bolometric flux peaks at redder wavelengths, increas- 
ing the flux in the NIR. In some cases one can also see the 
slight dip in the number of stars jus t brig htward of the 

(for example. 



TRGB that was discussed in Section 3.2.1 



see UGC4459, NGC4163, NGC3741 NTJG2976-DEEP. 
and IG2574-SGS for particularly obvious cases). 

To quantify the variation in the TRGB, we measure the 
TRGB magnitude in F160W using the edge-detection 
filter described in Mendez et al. (2002) a pplied to a 
Gauss ian-s moothed luminosit y function as in |Sakai et al.| 



cated techniques exist (e.g., Makarov et al. 2006 



(|1996|) and |Seth et al.| (|200 5|. Although more sophisti 

iMs 



Frayn 



& Gilmore 20031, the TRGB of our sample is typically 



well-populated and falls well above the photometric limit 
of the data, making our use of the widely used and cal- 
ibrated edge-detection technique adequate for an initial 
distance measurement. We use the identical procedure 
as was used to measure the TRGB for the optical data 



(iDalcanton et al.||2009] ), making the F8UW and NIR 
I'RGB directly comparable. Specifically, after extinc- 
tion correcting all magnitudes and colors, we restrict our 
analysis to stars on the RGB by iteratively fitting a line 
to stars that are less than 1 magnitude fainter than the 
estimated TRGB and that have colors consistent with 
potential RGB stars (0.6 < FllOW-FieOW < 1.1). We 
retain all stars that are within 2a of the fit to the RGB 
sequence towards the red, and within 1.5a to the blue; we 
use the more restrictive cut towards the blue to suppress 
the contribution from RHeB stars. The candidate RGB 
stars were then used to construct a Gaussian-smoothed 
luminosity function, which was then passed through an 
edge-detection filter. 

The final TRGB magnitude and uncertainty was mea- 
sured by executing 750 Monte Carlo bootstrap resam- 
pling trials. In each trial, additional Gaussian random 
errors were added to the stars' photometry, based on 
the magnitude of each star's photometric error. The 
TRGB for each trial was taken to be the magnitude corre- 
sponding to the peak of the edge-detection response filter 
within a 1 magnitude interval around the likely TRGB. 
We then fit the histogram of the returned TRGB magni- 
tudes with a Gaussian at totrgBj taking the mean and 
width of the Gaussian to be the magnitude of the TRGB 
and its uncertainty. These quantities are converted to ab- 
solute magnitudes using the distance moduli in Table [l] 
adopted from the TRGB analysis in FSUW. 

Once the magnitude of the TRGB has been measured, 
we characterize the color of the TRGB using stars within 
0.05 magnitudes fainter than the TRGB; in the few cases 
where there are fewer than 15 stars within the adopted 
magnitude range, we increase the range to 0.1 magni- 
tudes. We calculate the mean and standard deviation of 
these stars' colors using a biweight, iteratively clipping 
stars within 2a of the biweight mean, which visual in- 
spection suggests captures all the width of the TRGB 
with acceptable rejection of RHeB stars. 

Figure 18 shows the resulting absolute TRGB magni- 
tude in i^l60M^, as a function of the color of the TRGB. 
The data show the expected trend of brighter TRGB 
magnitudes for redder colors. This correlation can be 
well approximated as 



AfTRGB(i^l60W^) = -2.576(F110iy - F160W) - 3.496 

(1) 

with the absolute magnitudes showing a scatter of 0.050 
magnitudes around the mea n. 
Also shown on Figure 



18 are expectations for two sets 



of theoretical models. The red line shows the magnitude 
and color of the TRGB for 10 Gyr old Pado va isoc hrones 
from the WFC3 extension of |Girardi et al. (20081, span- 
ning a range of metallicities. I'he blue line shows a fit 
to the magnitude and color of the TRGB for a series of 
isochrones supplied by Aaron Dotter, spanning a range of 
ages (6-14 Gyr, shown as asteri with increasing point sizes 
for older TRGBs) and metalhcities (from [Fe/H]=-2.49 
(magenta) to [Fe/H]= —0.99 (green)). In general, the 
models show the same trend between TRGB color and 
magnitude seen in our data, suggesting that the bulk of 
the TRGB luminosity variation is driven by variations in 
metallicity. Similar behavior has been hinted at longer 
wavelengths for Spitzer IRAC observations of galaxies 
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within the Local Group ( Boyer et al. 2009 ) 



Although the slope of the ht to the Uotter isochrones is 
nearly identical to a weighted (in X and Y) least squares 
fit to the data (Eqn.fTl), the measured TRGBs are shifted 
to brighter magnitudes by a median of 0.099±G.027 mag- 
nitudes (or equivalently, to bluer colors by 0.043± 0.012), 
compared to the models. These offsets between the data 
and models are statistically significant, as judged by 
whether or not the mean residual in magnitude or color 
is significantly different from zero. 

We have considered a number of possible origins for 
the offset between data and models in the left hand 
plot of Figure [TSj First, as a broad check on the con- 
sistency of our measurements, on the right side of Fig- 
urejlSjwe plot the correlation between our measured NIR 
TRGB color {FllOW - FimW) and the optical-NIR 
color of the TRGB inferred from the difference between 
our measurements of mxRCB in F?>1AW and FIQOW . 
These two colors appear to be highly correlated, as ex- 
pected (Figure [T5| ) . However, they also show a small 
offset from the Dotter isochrones, such that the median 
observed optical-NIR color of the TRGB is 0.04 magni- 
tudes redder than the Dotter isochrones, with a semi- 
interquartile range of ±0.037. If the source of the offset 
in the MTRGB(i^l60VF) vs FllQW - Fl&QW relation 
were primarily a bias towards bluer measurements of the 
mean color of the TRGB, then the data should be offset 
to bluer colors in this diagram as well, rather than to 
redder FllOM^ — FIQQW colors as is observed. Thus, 
the offset between data and models does not appear to 
be due a mismeasurement of the TRGB color. 

Alternatively, if the source of the offset were primar- 
ily a bias towards brighter measurements of the FIQQW 
TRGB magnitude, then the data would be shifted to red- 
der values of to(F814VI^) — m{F\&QW) in this diagram. 
We do indeed see a redward shift between the data and 
the models for m(F814W^) — m(F160W^). However, the 
amplitude of this shift is a factor of 2.5 times smaller than 
needed to explain the offset observed in Mtrgb [FIQQW) 

vs FllOW -Fimw. 

The final explanation we consider is if the offset re- 
sults from a mismatch between the filter calibrations 
and the filter throughputs adopted by the models (such 
that measured F160W magnitudes are brighter than 
would be inferred for the models). If so, then both the 
m{FSUW) - m(F160W^) and FllQW - FimW values 
would be shifted to redder colors compared to the mod- 
els. However, this shift would move points along the 
mean relation, rather than perpendicular to it, and thus 
is unlikely to produce significant offsets. This possibility 
is consistent with the direction and the magnitudes of 



Fimw-H 



and 



0.2031-1-0.401 (J-i/-0.9)+0.3498 {J~H~0 

(2) 



J-H = 0.9418+0.841 (F110W-i^l60lV-1.0)-0.9053 {FllOW ~-F160W 

(3) 

which were derived by fitting the magnitudes and colors 
of the TRGB for colors redder than FUOW - F160W < 
1.15 in the Padova models; these transformations are 
good to < 0.001 magnitudes across this range. 

The globular cluster data tend to be shifted to brighter 
magnitudes and/or bluer colors than our measurements, 
and are even more offset from the models (in both the 
WFC3/IR and native J + H filter sets). The offsets are 
more pronounced for the Valenti sample, which appears 
to be biased somewhat brightwards compared to other 



data, includi n g two clusters in common with the Ivanov 
& Borissova (2002) sampl^^ and three bulge cl usters 



m common with the sample of Chun et al. (2010). We 
note, however, that measurements of the 'i'RGB can be 
difficult in globular clusters, which frequently have uncer- 
tain distances, sparsely populated RGBs at bright mag- 
nitudes, and ambiguous distinctions between bright RGB 
and faint AGB stars. Furthermore, the extinction correc- 
tions are typically quite large for the comparison globu- 
lar cluster sample; even with our restriction on the fore- 
ground extinction, the median extinction is larger than 
0.3 magnitudes in H for the Valenti sample. 

We have also considered that photometric biases may 
be responsible for the small offsets in our TRGB mag- 
nitudes and/or colors. For example, the Monte Carlo 
process used to evaluate our uncertainties artificially in- 
creases the photometric error (during randomization of 
magnitudes) and potentially biases totrgb by scattering 
stars preferentially above the tip. In practice, however, 
the effect of the added noise is negligible, since the photo- 
metric uncertainties are extremely small at ttitrgb ■ The 
effect of these biases have b een explored extensively in 



the offsets in both diagrams shown in Figure 18 



Madore & Freedman ( 19951, and at our typical signal- 
to-noise level of K ~ 50, we expect negligible magni- 
tude biases due to photometric errors. We likewise have 
considered the effect of crowding on our measurements, 
and find that it too is unlikely to explain the offset. As 
shown in Figure [Sj our photometric measurements are 
essentially unbiased in the median, with less than an 
0.005 magnitude shift at faint magntidues. Crowding 
does produce a slight tail to brighter magnitudes, due to 
undetected stars contaminating the flux of the resolved 
stars. However, at the typical magnitude of the TRGB, 
fewer than 3% of fake stars have magnitudes that would 
be brighter than expected from photometric uncertainty 
alone. In addition, we see no correlation between the de- 
gree of offset from the models and the number of stars 
in the field. 

We also have explored population differences as being 

The two clusters are NGC 6 441, which is . 07 ma g brighter in 
F160W and 0.08 mag bluer in the |Valenti et al.| | |2004l l sample, and 
NGC 6624, which is 0.12 m agii itudes brignter and Has an identical 
color in the Valenti et al. (2004[l sample. Both of these clusters have 
relatively high degrees ot toreground extinction, which increases the 
Including updates from http : //www. bo . astro . lt/~GC/ir_archlve/likelihood of differences between independent analyses. 



As a second approach to assessing the origin of the 
offset, we have included measurements of the TRGB from 
glob ular clusters in Figur e [T8l T hese measurements are 
from Ivanov fc Borissova (|2002 ) (s olid green diamonds) 
and the [Valenti et al.| ^'loU^ and [Valenti et al] ( [2007 1 
samples (open yellow diamond^^ restricting the sample 
to those with Ah < 0.5). Both oata sets were originally 
taken in J and H. We transformed these data to F160W 
and i^llOVF using the following relations 
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a source of the offset between the data and the models. 
Variations in the mean stellar age can affect the mag- 
nitude of the TRGB, while contamination from RHeB 
and AGB stars can change both the color and magnitude 
of the TRGB. However, we found no significant impact 
from these effects, based on the lack of correlation be- 
tween the the amplitude of the offset and the fraction of 
stars formed in different age ranges. 

We are left without any satisfying explanation for this 
small offset between the data and the models. At this 
point, our best guess is that it results from some com- 
bination of small uncertainties in the current WFC3/IR 
calibrations (although these improved significantly over 
the course of this program), in the filter throughputs 
adop ted by the mo dels, or the models themselves. In- 
deed, Cassisi (2010) Figure 2 shows significant systematic 
variations among models for the predicted bolomctric 
magnitudes of the TRGB. For now, any attempt to use 
F160W measurements for the TRGB should acknowl- 
edge that there may be an underlying 10% systematic 
uncertainty in the calibration. 

5. LUMINOSITY FUNCTIONS 

In the above analysis, we identify the TRGB as the 
point where there is a sharp peak in the edge-detection 
algorithm that also corresponds with the largest drop in 
the number of stars. However, the edge-detection algo- 
rithm frequently identifies other sharp transitions in the 
luminosity function that correspond to smaller changes 
in the absolute number of stars with magnitude. 

To explore these effects, in Figure[T9]we plot all galaxy 
luminosity functions, relative to the observed magnitude 
of the TRGB. The light lines are color coded accord- 
ing to each galaxy's rank when sorted by the fraction of 
star formation the galaxy has experienced in the most 
recent Gyr. The heavy lines show the average luminos- 
ity function when the galaxies are sorted by rank into 
4 bins of recent star formation, with equal numbers of 
galaxies per bin; red, magenta, blue, and black lines go 
from lowest to highest fraction of recent star formation 
((/o-iGyr) = 0.001,0.04,0.07,0.09, respectively). AU 
luminosity functions have been normalized to have the 
same total number of stars in the bins between 0.25 and 
0.75 magnitudes fainter than the TRGB. We do not ana- 
lyze the luminosity function at fainter magnitudes, where 
the varying depths, photometric uncertainties, and dis- 
tances among the galaxy population make the averag- 
ing procedure invalid. These luminosity functions in- 
clude only stars that lie within the extrapolated slope 
and width of the RGB, and may miss the reddest, most 
luminous RHeB stars. 



The average luminosity functions in Figure 19 show 
clear systematic deviations that correlate with the frac- 
tion of star formation in the most recent gigayear. These 
deviations are most obvious brighter than the TRGB, 
where elevated recent star formation produces a larger 
population of RHeB and luminous AGB stars, extend 
ing to brig hter magnitudes. As not ed previously in Sec 



tion |3.2.2 and as quantified in , Melbourne et al. (2011), 
this population of luminous stars can lead to signin- 
cant reductions in the NIR mass-to-light ratio, on much 
shorter timescales than discussed for AGB stars alone. 

There are also subtle features close to the TRGB it- 
self, some of which correlate with the recent SFR. In 



particular, the amplitude of the discontinuity at the 
TRGB is less pronounced when recent star formation is 
more prominent, most likely due to contamination from 
younger stellar populations all along the RGB sequence 
and its extension to brighter magnitudes. This contam- 
ination will tend to reduce the reliability of the NIR 
TRGB distance measurements for actively star-forming 
galaxies. 

6. METALLICITY AND THE STRUCTURE OF THE RGB 

For two decades it has been known that the slope of the 
NIR RGB correla tes strongly with metallicity for globu- 
lar clusters (e.g., 'Davidgc e t aL] 1992[ Cohen & Sleeper 
19951 iKuchinski et al. 1995 ; iJj'erraro et al.||'2000|). WEIe 



this correlation allows the slope of the INIR RGB to be 
used as a reddening-free diagnosti c of metallicity for gl ob 



ular clusters (for an example, see Ferraro et aT] 2006 1 , it 
is not clear if the relationships used for uni-aged globular 
clusters can be routinely applied to the complex stellar 
populations found in galaxies, where age variations can 
also affect the structure of the RGB, and where observa- 
tions typically produce lower quality photometry due to 
unavoidable crowding errors. 

We now evaluate whether the NIR RGB slope can be 
used effectively as a metallicity indicator for complex 
stellar populations. We first characterize the RGB using 
the measured median color as a function of magnitude 
(red lines from Figure 17 see Section |4|. We presume 
that the observed variation in color is dominated by vari- 
ations in metallicity, based on the greater sensitivity of 
RGB color to metallicity compared to age (upper right 
of Figure 10). Thus, while we cannot directly compare 
the RGB slope to metallicity, we can use the color of the 
TRGB as a proxy. 

6.1. TRGB Color and the Luminosity- Metallicity 
Relationship 

To demonstrate the connection between metallicity 
and NIR color of the TRGB, in Figure [20| we plot 
the color of the TRGB as a function of each galaxy's 
extinction-corrected magnitude (upper left) and NIR 
luminosity (upper right; 3.6/i luminosity taken from 
(Dale et al. 20091, assuming that the Sun has a flux 
of 14.71 Jy at a distance of 10 pc). The data show a 
clear relationship between galaxy luminosity and TRGB 
color,such that FllOW ~ F160W = 0.605 - 0.0181 A/b^ 
and FUOW - F160W = 0.431 + 0.0518 log^o Lg.g^, with 
an rms of 0.046 and 0.039, respectively. Beca use of the 
well-k nown mass-metallicity relationship (e.g., Lee et al. 
20061 and references therein), one would expect higher 



metallicities in more massive, luminous galaxies, which 
should then manifest itself as mor e luminous galaxies 
having redder TRGBs. Figure 20 shows this expected 
behavior, suggesting that the TRGB color does indeed 
correlate with metallicity. 

Further support for the connection between TRGB 
color and metallicity comes from the lower left panel, 
where we plot the NIR TRGB color as a function of 
the metallicity inferred from the oldest age bin of the 
star formation histor y (derived from the opti cal photom- 
etry using MATCH; jMelbourne et al.||2011| . We see a 
tight correlation between JNIR I'RGB color and metal- 
licity, following the relationship FllOW — i^l60M^ — 
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1.093 + 0.192 [Fe/H] with an rms of 0.050. Note, how- 
ever, that one expects this correlation to be tight, given 
that the metallicity sensitivity of the CMD fits conies 
largely from the color of the TRGB. 

6.2. RGB Slope as a Metallicity Indicator 

Having verified the relationship between TRGB color 
and metallicity, we now investigate the correlation be- 
tween RGB slope and the TRGB color. On the left side of 
Figure [21] we plot linear fits to the RGB, restricted to two 
magnitudes fainter than the TRGEr^ The plot shows the 
expected variation of the magnitude of the TRGB with 
color, such that redder RGBs extend to brighter magni- 
tudes (Figure 18). 



On the right side of Figure [21] we plot the slope of the 
linear fit compared to the color of the fit at the magni- 
tude of the TRGB. In general, there are no systematic 
statistically-significant trends across the whole sample. 
There is a ^25% dispersion in the mean slope for colors 
between 0.8 < FllOW - F160W < 1, with no apparent 
correlation. This suggests that the observed RGB slope 
is unlikely to be a useful metallicity indicator in the pres- 
ence of the dispersion in age, metallicity, and photomet- 
ric errors within individual extragalactic systems. 

The only possible manifestations of the known corre- 
lation between NIR RGB slope and metallicity can be 
seen at the extremes. Of the four galaxies with the 
most vertical RGBs, three also have the bluest colors. 
Two of these three (UA292 and Sc22) do indeed to have 
very low metallicity populations, based upon the lack 
of curvature in the optical RGB. In contrast, the re- 
maining galaxy (KDG73) does not have a particularly 
old globular-like stellar population, nor does the optical 
CMD suggest a particularly low metallicity (based upon 
the observed curvature in the optical RGB). Instead, this 
galaxy appears to have a substantial younger AGB pop- 
ulation which may be pulling the RG B to bluer colors 
and altering its morphology (see Figure 10 1. 

The other extreme outliers are HoII ana IC2574, both 
of which have particularly shallow RGB slopes. However, 
both of these galaxies also have dramatic RHeB popula- 
tions, which are likely shifting the base of the RGB fit to 
bluer colors, artificially flattening the RGB slope. 

In summary, there appears to be little chance that the 
slope of the NIR RGB can be used to assess the metallic- 
ity of extragalactic systems with current data quality; the 
existing trends appear to be weak at best, and outliers 
due to contributions from AGB and RHeB stars are not 
uncommon. One may have better success by combining 
the NIR with the optical to produce an optical-NIR CMD 
with a wide color baseline (i.e., similar to the widely used 
V — K color) , or with using a more sophisticated method 
to characterize the RGB slope. 

7. CONCLUSIONS 

This work represents the first step in characterizing the 
NIR properties of complex stellar populations for galax- 
ies spanning a wide range of metallicity and star forma- 
tion history. We have presented CMDs in one of the 

The RGB slope in globular clusters is typically measured down 
to near the horizontal branch (e.g., Kuchinski et al. 1995). How- 
ever, in extragalactic systems one typically cannot resolve more 
than 2-3 magnitudes of the RGB, due to crowding. 



most efficient and commonly used WFC3/IR filter pairs, 
allowing one to assess the contribution of individual stars 
to the integrated light in these bandpasses. 

In all cases, the CMDs show a dominant population 
of RGB stars, as expected. They also show clear se- 
quences of AGB stars, with only a modest contribution 
of extremely red AGB stars. These features are charac- 
teristic of old and intermediate age star formation. In 
galaxies with more recent star formation, we also find 
a dramatic sequence of red core Helium burning stars, 
extending to much brighter magnitudes than the AGB 
sequence. These stars are associated with recent star 
formation (20-500 Myr), and in some cases are a major 
contributor to the NIR flux. They also have colors that 
are similar to the underlying AGB and RGB population, 
making their presence difficult to diagnose from broad- 
band colors alone. These RHeB stars therefore present a 
significant uncertainty in adopting a NIR mass-to-light 
ratio when analyzing unresolved galaxies or interpreting 
their NIR colors. We qu antify the contribution of both 
RHeB and AGB stars in [Melbourne et al.| ( [20TT] ). 

We have used the observed CMDs to empirically cal- 
ibrate the NIR TRGB magnitude as a function of the 
metallicity-sensitive FllOW — F160W color, allowing 
these standard HST filters to be used as distance in- 
dicators in the WFC3/IR and JWST era. We find that 
there is a strong correlation between color and the abso- 
lute magnitude of the TRGB, as expected from models 
and previous observations of globular clusters. However, 
we find that the relationship is offset by ~0. 05-0.1 mag- 
nitudes from the prediction of current isochrone models. 
Analysis of the origin of this offset points to a residual 
uncertainty in either the WFC3/IR zero points or the 
theoretical models. 

We also explore the structure of the RGB and AGB se- 
quences, highlighting possible variations with the age and 
metallicity of the underlying stellar population. We see 
a clear age-dependent variation in the luminosity func- 
tion of red stars, due to increasing numbers of luminous 
AGB stars with increasing intermediate age star forma- 
tion. We are unable to detect any metallicity dependence 
in the slope of the RGB, however, which prevents the 
slope of the NIR RGB from being used as a reddening- 
free metallicity indicator. 

We also present an analysis of scattered light in the 
WFC3/IR detector. We have found a rapidly- varying 
scattered light component that affects images taken at 
low angles to a bright Earth limb. We discuss how this 
scattered light may affect photometry, particularly when 
fitting a series of non-destructive reads with a linear func- 
tion. 

In a series of subsequent papers, we will be using the 
data presented here to refine theoretical models of RHeB 
and AGB stars, to assess contributions of carbon stars, 
to empirically calibrate the fractional luminosity due to 
AGB and RHeB stars, and to improve star formation 
histories at intermediate ages. 
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TABLE 1 
Sample Galaxies 



Galaxy 


Alt. 


RA 


Dec 


Diam. 




Ay 


m - M 


T 




Group 




Names 


(J2000) 


(J2000) 


(') 














DD053 


U4459 


08:34:06.5 


66:10:45 


1.6 


14.55 


0.118 


27.79 


10 


25 


M81 


DD078 




10:26:27.9 


67:39:24 


2.0 


15.8 


0.066 


28.18 


-3 




M81 


DD082 


U5692 


10:30:35.0 


70:37:10 


3.4 


13.57 


0.133 


27.90 


9 




M81 


Hoi 


U5139,DD063 


09:40:28.2 


71:11:11 


3.6 


13.64 


0.153 


27.95 


10 


29 


M81 


HoII 


U4305 


08:19:05.9 


70:42:51 


7.9 


11.09 


0.098 


27.65 


10 


66 


M81 
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10:21:25.2 


71:06:58 


1.5 


16.5 


0.359 


27.91 


10 


13 


M81 


12574 


U5666,DD081 


10:28:22.4 


68:24:58 


13.2 


10.84 


0.112 


27.90 


9 


115 


M81 


KDG2 


E540-030.KK9 


00:49:21.1 


-18:04:28 


1.2 


16.37 


0.072 


27.61 


-1 




Scl 


KDG63 


U5428,Df)071 


10:05:07.3 


66:33:18 


1.7 


16.01 


0.303 


27.74 


-3 


19 


M81 


KDG73 




10:52:55.3 


69:32:45 


0.6 


17.09 


0.056 


28.03 


10 


18 


M81 


KKH37 




06:47:45.8 


80:07:26 


1.2 


16.4 


0.231 


27.56 


10 


20 




M81 


N3031,U5318 


09:55:33.5 


69:04:00 


26.9 


7.69 


0.249 


27.77 


3 


422 


M81 


N300 




00:54:53.5 


-37:40:57 


21.9 


8.95 


0.039 


26.5 


7 


149 


14-(-13 


N404 


U718 


01:09:26.9 


35:43:03 


2.5 


11.21 


0.181 


27.42 


-1 


78 




N2403 


U3918 


07:36:54.4 


65:35:58 


21.9 


8.82 


0.124 


27.5 


6 


231 


M81 


N2976 


U5221 


09:47:15.6 


67:54:49 
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11.01 


0.224 


27.76 


5 


97 


M81 


N3077 


U5398 


10:03:21.0 


68:44:02 


5.4 


10.46 


0.208 


27.92 


10 


65 


M81 


N3741 


U6572 


11:36:06.4 


45:17:07 
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14.38 


0.077 


27.55 


10 


81 


14-(-07 


N4163 


U7199 


12:12:08.9 


36:10:10 


1.9 
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0.062 


27.29 


10 


18 


li+07 


N7793 




23:57:49.4 


-32:35:24 


9.3 


9.70 


0.060 


27.96 


7 
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Sc22 


Sc-dEl 


00:23:51.7 
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17.73 


0.046 


28.11 


10 




Scl 
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IZw60 


13:30:44.4 


54:54:36 


1.7 


14.12 


0.047 


27.06 


10 


49 


14-1-07 


UA292 


CVnI-dwA 


12:38:40.0 


32:46:00 
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16.10 


0.048 


27.79 


10 


27 





Note. — Distances, -Bt, and T-type taken fr om CNG, with updat es f rom l^ alcaiiton et al .|2009[ l and jKarachentsev| 

|et al.|2003l l for NGC 7793; Group membership from |Karachentsev| j2005[ or |'ltilly et aL| | |2UUe^ ; Ay from mfc^ 



TABLE 2 
Observations & Photometry 



Catalog 


Target 




Obs.Date 




N stars 






50% Completeness 


Optical 


Optical 


Name 


Name 














(#/sqr- 


■arcsec) 


(FllOW) 
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PropID 


Filters 


DD053 


UGC4459 


2010- 


■04-23 


11: 


:15: 


:02 


9366 


2.24 


0.07 


26.02 


25.02 


GO-10605 


F555W,F814W 
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2010- 
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2010- 
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2010- 
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27898 
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2010- 
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2010- 
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F606W,F814W 
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SN-NGC2403-PR 


2010- 
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5.77 
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F475W,F606W,F814W 


N2976 
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2010- 
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15392 


4.05 


0.25 


25.84 


24.92 


GO-10915 


F475W,F606W,F814W 


N3077 


NGC3077-PHOENIX 


2010- 


■02-21 


22: 


:49: 


:04 


8813 


3.14 


0.20 


26.00 


25.14 


GO-9381 


F435W,F555W,F814W 


N3741 


NGC3741 


2009- 


■11-07 


01: 


:33: 


:13 


6819 


4.44 


0.00 


25.98 


24.86 


GO-10915 


F475W,F814W 


N4163 


NGC4163 


2010- 


■03-23 


17: 


:39: 


:56 


19105 


4.77 


0.07 


25.90 


24.71 


GO-10915 


F475W,F606W,F814W 


N7793 


NGC7793-HALO-6 


2010- 


■06-14 


19: 


:11: 


:39 


6578 


2.14 


0.02 


26.11 


25.07 


GO-10523 


F606W,F814W 


Sc22 


SCL-DEl 


2009- 


■09-08 


00: 


:45: 


:13 


2440 


3.17 


0.02 


26.19 


25.09 


GO-10503 


F606W,F814W 


U8508 


UGC8508 


2009- 


■10-14 


19: 


:41: 


:41 


12664 


4.38 


0.03 


25.97 


24.83 


GO-10915 


F475W,F814W 


UA292 


UGCA292 


2010- 


■05-18 


12: 


:38: 


:27 


2373 


2.53 


0.01 


26.19 


25.14 


GO-10915 


F475W,F814W 



Note. — All star counts and surface densities are based on the numbers of stars in the *.gst catalogs; the * . st catalogs typically contain twice as many 
stars, though with larger photometric uncertainties. The maximum and minimum stellar surface densities were calculated by calculating the density of sources 
in a 10x10 grid on the image. 
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TABLE 3 
TRGB Measurements 



Catalog 


Target 


m — M 


Ay 


stars 


Mean 




rriTRGB 




Name 


Name 


(F814VK) 






(TRGB fit) 


Color 


(raw) 


{F160W) 


(F160iy) 




DD071 


27.74 


0, 


.303 


1024 


0.853 


22.142 


22. U80 ± U.064 


— 5.660 ± 0.064 


DD078 


DDU78 


27.82 


0, 


.066 


1438 


0.901 


22.045 


OO 00"l 1 f\ /~\o^ 

22.031 ± 0.027 


—5.789 ± 0.027 




DDU82 


27.90 


0, 


.133 


5268 


0.909 


22.089 




CT ooo 1 /~i noo 

— o.ooo ± U.U25 




-tbUo4U-UoU 


27.61 


0, 


.072 


481 


0.829 


22.105 


Zl.Um) ± U.U4U 


— O.OzU ± U.U4U 


HS117 


HS117 


27.91 


0, 


.359 


480 


0.836 


22.340 


22.2bb ± (J.Uzo 


—5.644 ± 0.025 


12574 


i02574-bGb 


27.90 


0, 


.112 


3967 


0.859 


22.232 


oo or»ri 1 o r\i r\ 

22.2\yd ± (J.U19 


— o.D9i ± U.U19 


KJJCj / o 




28.03 


0, 


.056 


273 


0.771 


22.464 


ZZ.40Z ih U.U40 


— 0.0 r o ± U.U4o 




KKH37 


27.56 


0, 


.231 


659 


0.848 


21.957 


21.910 ± 0.026 


cr c n 1 n no^ 

—5.650 ± 0.026 


M81 


M81-UrjrjF 


27.77 


0, 


.249 


594 


0.986 


21.744 


21.o9o ± (J.U43 


— D.U77 ± (J.U4^ 


N300 


IN CjUUoUU- W iUrji 


26.50 


0, 


.039 


1410 


0.972 


20.561 


ZU.OOO ± U.UzD 


— 0.94 r ± U.Uzo 


MO /I r\Q 

JNz4Uo 


IN GC24Uo-DrvrjF 


27.50 


0, 


.124 


751 


0.927 


21.701 


oi cjTcr 1 or\o£; 

21.675 ± 0.036 


—5.825 ± 0.036 


JNz4Uo 


1nCj024Uo-HALU-0 


27.50 


0, 


.124 


356 


0.901 


21.635 


21.6U9 ± 0.016 


— 5.891 ± 0.016 


MOriT/^ 

JN^97d 


IN CjOzy 7o-Drjrjl^ 


27.76 


0, 


.224 


2495 


0.952 


21.851 


21.8U5 ± 0.034 


— 5.955 ± 0.034 






27.92 


0, 


.208 


945 


1.000 


21.931 


21.888 ± 0.020 


— D.U32 ± 0.025 


N3741 


NGC3741 


27.55 


0, 


.077 


801 


0.814 


21.974 


21.958 ±0.023 


-5.592 ±0.023 


N404 


NGC404 


27.42 


0, 


.181 


2691 


0.967 


21.498 


21.461 ± 0.021 


-5.959 ±0.021 


N4163 


NGC4163 


27.29 


0, 


.062 


2594 


0.875 


21.571 


21.558 ±0.014 


-5.732 ±0.014 


N7793 


NGC7793-HALO-6 


27.96 


0, 


.060 


648 


0.949 


21.930 


21.917 ± 0.020 


-6.043 ±0.020 


Sc22 


SCL-DEl 


28.11 


0, 


.046 


377 


0.793 


22.570 


22.560 ± 0.037 


-5.550 ± 0.037 


HoII 


UGC4305-1 


27.65 


0, 


.098 


1845 


0.850 


21.937 


21.917 ± 0.071 


-5.733 ±0.071 


HoII 


UGC4305-2 


27.65 


0, 


.098 


2258 


0.835 


21.947 


21.927 ± 0.025 


-5.723 ±0.025 


DD053 


UGC4459 


27.79 


0, 


.118 


932 


0.805 


22.142 


22.118 ±0.031 


-5.672 ±0.031 


Hoi 


UGC5139 


27.95 


0, 


.153 


1742 


0.826 


22.351 


22.320 ± 0.029 


-5.630 ± 0.029 


U8508 


UGC8508 


27.06 


0, 


.047 


1110 


0.847 


21.412 


21.402 ±0.024 


-5.658 ±0.024 


UA292 


UGCA292 


27.79 


0, 


.048 


184 


0.722 


22.459 


22.450 ± 0.022 


-5.340 ± 0.022 



Note. — Absolute magnitudes arc determined using the distance modulus from Table |^ which were originally determined from 
the FS14:W TRGB. Mean colors are for the stars used to measure the TRGB, which arc not necessarily all RGB stars, and include 
only stars within 0.05"^ of the TRGB. Av values are as reported by IRSA for coordinates in Table 1. Extinction corrections from 
Av to the observed filters are adopted from models of |Girardi et ah] ( |2008j l, as described in the text. Listed uncertainties are 
dominated by photometric uncertainties and by stochasticity in the number of stars near the tip; systematic uncertainties (due to 
uncertainties in the assumed TRGB absolute magnitudes and extinction) are likely to be much larger, but are not included in the 
listed uncertainties. SN-NGC2403-PR is not included due to large crowding errors. 





Fig. 1. — Position of the WFC3/IR field of view (red), overlayed on an optical image from the Digitized Sky Survey. Blue region shows 
the area covered by optical HST data. (Target names from upper left to lower right: [a] UGC4459; [b] DD078; [c] DD082; [d] UGC5139; 
[e] UGC4305-1; [f] UGC4305-2; [g] HS117; [h] IC2574-SGS; [i] ESO540-030;). 
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Fig. [T] CONTINUED.— Position of the WFC3/IR field of view (red), overlayed on an optical image from the Digitized Sky Survey. Blue 
region shows the area covered by optical HST data. (Target names from upper left to lower right: [j] DD071; [k] KDG73; [1] KKH37; [m] 
M81-DEEP; [n] NGC0300-WIDE1; [o] NGC404; [p] NGC2403-DEEP; [q] NGC2403-HALO-6; [r] SN-NGC2403-PR;). 





Fig. 2.— (Left) False color FUOW + FWOW image of the WFC3/IR field for the target UGC4459 within the galaxy DD053 [Top] and 
within the galaxy DD078 [Bottom]. (Right) Color- magnitude diagrams generated for a grid of subregions, so that the upper left CMD 
corresponds to the upper left of the adjacent image. 




Fig. [2] CONTINUED.— (Left) False color FllOW + F160W image of the WFC3/IR field for within the galaxy DD082 [Top] and the target 
UGC5139 within the galaxy Hoi [Bottom]. (Right) Color-magnitude diagrams generated for a grid of subregions, so that the upper left 
CMD corresponds to the upper left of the adjacent image. 
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Fig. \2\ CONTINUED.— (Left) False color FUOW + F160W image of the WFC3/IR field for the target UGC4305-1 within the galaxy HoII 
[Top] and the target UGC4305-2 within the galaxy HoII [Bottom]. (Right) Color-magnitude diagrams generated for a grid of subregions, 
so that the upper left CMD corresponds to the upper left of the adjacent image. 
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Fig. [1 CONTINUED.— (Left) False color FllOW + F160W image of the WFC3/IR field for within the galaxy HS117 [Top] and the target 
IC2574-SGS within the galaxy 12574 [Bottom]. (Right) Color-magnitude diagrams generated for a grid of subregions, so that the upper 
left CMD corresponds to the upper left of the adjacent image. 
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Fig. [2] CONTINUED.— (Left) False color FllOW + FWOW image of the WFC3/IR field for the target ESO540-030 within the galaxy KDG2 
[Top] and the target DD071 within the galaxy KDG63 [Bottom]. (Right) Color-magnitude diagrams generated for a grid of subregions, so 
that the upper left CMD corresponds to the upper left of the adjacent image. 
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Fig. [2] CONTINUED.— (Left) False color FllOiy + FieOW image of the WFC3/IR field for within the galaxy KDG73 [Top] and within the 
galaxy KKH37 [Bottom]. (Right) Color-magnitude diagrams generated for a grid of subregions, so that the upper left CMD corresponds 
to the upper left of the adjacent image. 
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Fig. [2] CONTINUED.— (Left) False color i^'llOW' + F160VF image of the WFC3/IR field for the target M81-DEEP within the galaxy 
M81 [Top] and the target NGC0300-WIDE1 within the galaxy NSOO [Bottom]. (Right) Color-magnitude diagrams generated for a grid of 
subregions, so that the upper left CMD corresponds to the upper left of the adjacent image. 
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Fig. [2] CONTINUED.— (Left) False color FllOW + F160W image of the WFC3/IR field for the target NGC404 within the galaxy N404 
[Top] and the target NGC2403-DEEP within the galaxy N2403 [Bottom]. (Right) Color-magnitude diagrams generated for a grid of 
subregions, so that the upper left CMD corresponds to the upper left of the adjacent image. 
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Fig. [1 CONTINUED. (Left) False color FllOW + F160H'' image of the WFC3/IR field for the target NGC2403-HALO-6 within the 
galaxy N2403 [Top] and the target SN-NGC2403-PR within the galaxy N2403 [Bottom]. (Right) Color-magnitude diagrams generated for 
a grid of subrcgions, so that the upper left CMD corresponds to the upper left of the adjacent image. 
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Fig. [2] CONTINUED.— (Left) False color FllOW + FWOW image of the WFC3/IR field for the target NGC2976-DEEP within the galaxy 
N2976 [Top] and the target NGC3077-PHOENIX within the galaxy N3077 [Bottom]. (Right) Color-magnitude diagrams generated for a 
grid of subregions, so that the upper left CMD corresponds to the upper left of the adjacent image. 
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Fig. [3 CONTINUED. (Left) False color FllOW + F160W image of the WFC3/IR field for the target NGC3741 within the galaxy N3741 
[Top] and the target NGC4163 within the galaxy N4163 [Bottom]. (Right) Color-magnitude diagrams generated for a grid of subregions, 
so that the upper left CMD corresponds to the upper left of the adjacent image. 
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Fig. [2] CONTINUED.— (Left) False color FllOW^ -I- FimW image of the WFC3/IR field for the target NGC7793-HALO-6 within the 
galaxy N7793 [Top] and the target SCL-DEl within the galaxy Sc22 [Bottom] . (Right) Color-magnitude diagrams generated for a grid of 
subregions, so that the upper left CMD corresponds to the upper left of the adjacent image. 
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Fig. [3 CONTINUED.— (Left) False color FUOW + F160VK image of the WFC3/IR field for the target UGC8508 within the galaxy U8508 
[Top] and the target UGCA292 within the galaxy UA292 [Bottom]. (Right) Color-magnitude diagrams generated for a grid of subregions, 
so that the upper left CMD corresponds to the upper left of the adjacent image. 
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Fig. 3. — Comparison between a false-color image in the optical (left: i<'435l'V' + ^555^' + F814W) and the NIR (right: F110W + F160W) 
in a star forming region of 12574. Nebulosity associated with line emission from ionized [Slll] and netural He is visible near Hll regions in 
the NIR image. The circle indicates a compact star formation region dominated by luminous red core-Helium burning (RHeB) stars. This 
SF region is old enough that the photoionized Hll has recombined (> lOMyrs), and therefore lacks the diffuse e mission seen in the younger 
neighboring HH regions. However, it also hosts all of the most NIR luminous sources (see discussion in Scction|3.2.2[l. 




Fig. 4. — DOLPHOT's estimate of the local sky for a sequence of exposures taken during a single orbit visit for UGC4305-1, based 
on DOLPHOT's PSF fitting. The first exposure in the orbit is in the upper left, and the final exposure is in the lower right. Exposures 
alternate between _F110VF and F160W, and all images for a single filter have been displayed with a common scale and stretch. There is 
an overall gradient from top to bottom in all images, due to a larger fraction of unresolved sources in the upper half of the image. There 
is also a time-dependent feature on the left hand side of the image, due to scattered light. We show in Figure [s] that this feature appears 
when the telescope is pointed to within 30° of a bright earth limb; the first two exposures in the orbit had average limb angles of 22° and 
26°, respectively, while the third through sixth exposures have limb angles of 31°, 34°, 36°, and 37°, respectively. Even after the bright 
limb feature disappears, there are still temporal variations in the overall sky brightness level. 
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Fig. 5. — (Left) Sky level (counts per pixel) in the upper right quadrant of all FllOW and Fl&OW exposures (top and bottom, 
respectively), as a function of the average angle between the spacecraft's VI axis and the earth's limb during the exposure; open circles 
indicate that the limb was dark at the time of observation. Successive exposures in a single orbit are connected by lines, but intermediate 
non-destructive reads are not shown. In FllOW , the mean sky level is clearly correlated with proximity to the earth's limb, when the limb 
is bright. (Right) The amplitude of the scattered light feature seen in Figure[4](as a fraction of the sky brightness), as a function of the 
average angle to the earth's limb during each observation. The amplitude of tne scattered light feature is taken to be the excess counts 
in the affected region, compared to the variation in the mean sky level seen in the unaffected upper right quadrant. It is calculated by 
first calculating a gradient across the image using the difference between the amplitude of the sky in the affected region and the unaffected 
region in the upper right of the image. The smallest gradient for the exposures in a given orbit is then taken to be the expected true spatial 
variation in the sky background, in the absence of scattered light. The amplitude of the scattered light feature is then scaled to the level 
of the sky in the upper right quadrant of each image, giving the approximate amplitude of the feature as a percent of sky. The scattered 
light feature turns on abruptly at earth limb angles less than 30°. 
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Fig. 6. — Photometric uncertainties as a function of magnitude in the FllOW (top) and F160W filters, for a representative field with 
moderate crowding (UGC8508), using the more complete *.st catalog. The vertical lines indicate the average 50% completeness level 
determined from artificial star tests (Table[2|. Photometric uncertainties are those reported by DOLPHOT, and do not fully capture the 
uncertainties due to crowding (see Figure 
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Fig. 7. — Comparison of the color magnitude diagrams of the original stellar *.st photometry (loft) and the high quality "cleaned" 
*.gst photometry (right) for target IC2574-SGS. The original photometry has higher completeness, but at the expense of increased errors 
due to crowding, producing less sharp features in the CMD. The *.st photometry includes stars that have a signal-to-noise of greater than 
4 in only one filter, making the colors unreliable at faint magnitudes; the *.gst photometry requires high signal-to-noise in both filters. 
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Fig. 8. — Distribution of differences between the true and the recovered F\&SW magnitudes for artificial stars inserted into images 
of IC2574-SGS (left; highly crowded) and UGC 4459 (right; uncrowded). Distributions are calculated in bins 1 magnitude wide, with 
increasing line thicknesses indicating fainter magnitudes; the heaviest line includes stars with 23 < mpieow < 24. Crowding biases the 
photometry very slightly (< 0.004 mag) towards fainter observed magnitudes at the median, with larger biases at fainter magnitudes. 
However, these biases are much smaller than the photometric errors (<0.05a-m for 1C2574-SGS) in all magnitude ranges. The uncertainties 
in the magnitude of individual stars is frequently dominated by crowding, however (see Figure^for typical photometric errors as a function 
of magnitude.) 
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Fig. 9. — Color magnitude diagrams of the WFC3/IR (upper left) and optical (upper right) for the target UGC4459 within galaxy 
DD053. Lower panels show the star formation history derived from the optical data, for both the differential (left, with horizontal dotted 
line indicating the past average SFR) and cumulative (right) star formation histories. The cumulative star formation history is calculated 
from the present back to 14 Gyrs. tFncertainties in the lower two panels are the 68% confidence intervals, calculated from Monte Carlo 
tests including random and systematic uncertainties. Optical CMDs are restricted to the area covered by the WFC3 FOV. 
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Fig. Is] CONTINUED. — Color magnitude diagrams of the WFC3/IR (upper left) and optical (upper right) for DD078. Lower panels show 
the star formation history derived from the optical data, for both the differential (left, with horizontal dotted line indicating the past 
average SFR) and cumulative (right) star formation histories. The cumulative star formation history is calculated from the present back 
to 14Gyrs. Uncertainties in the lower two panels are the 68% confidence intervals, calculated from Monte Carlo tests including random 
and systematic uncertainties. Optical CMDs are restricted to the area covered by the WFC3 FOV. 
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Fig. m CONTINUED. — Color magnitude diagrams of the WFC3/IR (upper left) and optical (upper right) for DD082. Lower panels show 
the star formation history derived from the optical data, for both the differential (left, with horizontal dotted line indicating the past 
average SFR) and cumulative (right) star formation histories. The cumulative star formation history is calculated from the present back 
to 14Gyrs. Uncertainties in the lower two panels are the 68% confidence intervals, calculated from Monte Carlo tests including random 
and systematic uncertainties. Optical CMDs are restricted to the area covered by the WFC3 FOV. 
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Fig. [H CONTINUED. — Color magnitude diagrams of the WFC3/IR (upper left) and optical (upper right) for the target UGC5139 within 
galaxy Hoi. Lower panels show the star formation history derived from the optical data, for both the differential (left, with horizontal 
dotted line indicating the past average SFR) and cumulative (right) star formation histories. The cumulative star formation history is 
calculated from the present back to 14Gyrs. Uncertainties in the lower two panels are the 68% confidence intervals, calculated from Monte 
Carlo tests including random and systematic uncertainties. Optical CMDs are restricted to the area covered by the WFC3 FOV. 
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Fig. [9] CONTINUED. — Color magnitude diagrams of the WFC3/IR (upper left) and optical (upper right) for the target UGC4305-1 within 
galaxy HoII. Lower panels show the star formation history derived from the optical data, for both the differential (left, with horizontal 
dotted line indicating the past average SFR) and cumulative (right) star formation histories. The cumulative star formation history is 
calculated from the present back to 14Gyrs. Uncertainties in the lower two panels are the 68% confidence intervals, calculated from Monte 
Carlo tests including random and systematic uncertainties. Optical CMDs are restricted to the area covered by the WFC3 FOV. 
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Fig. [9] CONTINUED. — Color magnitude diagrams of the WFC3/IR (upper left) and optical (upper right) for the target UGC4305-2 within 
galaxy HoII. Lower panels show the star formation history derived from the optical data, for both the differential (left, with horizontal 
dotted line indicating the past average SFR) and cumulative (right) star formation histories. The cumulative star formation history is 
calculated from the present back to 14Gyrs. Uncertainties in the lower two panels are the 68% confidence intervals, calculated from Monte 
Carlo tests including random and systematic uncertainties. Optical CMDs are restricted to the area covered by the WFC3 FOV. 





10.00 1.00 0.10 0.01 Age(Gyr) 
Age (Gyr) 

Fig. [9I CONTINUED. — Color magnitude diagrams of the WFC3/IR (upper left) and optical (upper right) for HS117. Lower panels show 
the star formation history derived from the optical data, for both the differential (left, with horizontal dotted line indicating the past 
average SFR) and cumulative (right) star formation histories. The cumulative star formation history is calculated from the present back 
to 14Gyrs. Uncertainties in the lower two panels are the 68% confidence intervals, calculated from Monte Carlo tests including random 
and systematic uncertainties. Optical CMDs are restricted to the area covered by the WFC3 FOV. 
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Fig. [U CONTINUED. — Color magnitude diagrams of the WFC3/IR (upper left) and optical (upper right) for the target IC2574-SGS within 
galaxy 12574. Lower panels show the star formation history derived from the optical data, for both the differential (left, with horizontal 
dotted line indicating the past average SFR) and cumulative (right) star formation histories. The cumulative star formation history is 
calculated from the present back to 14Gyrs. Uncertainties in the lower two panels are the 68% confidence intervals, calculated from Monte 
Carlo tests including random and systematic uncertainties. Optical CMDs are restricted to the area covered by the WFC3 FOV. 
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Fig. [U CONTINUED. — Color magnitude diagrams of the WFC3/IR (upper left) and optical (upper right) for the target ESO540-030 within 
galaxy KDG2. Lower panels show the star formation history derived from the optical data, for both the differential (left, with horizontal 
dotted line indicating the past average SFR) and cumulative (right) star formation histories. The cumulative star formation history is 
calculated from the present back to 14Gyrs. Uncertainties in the lower two panels are the 68% confidence intervals, calculated from Monte 
Carlo tests including random and systematic uncertainties. Optical CMDs are restricted to the area covered by the WFC3 FOV. 
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Fig. [U CONTINUED. — Color magnitude diagrams of the WFC3/IR (upper left) and optical (upper right) for the target DD071 within 
galaxy KDG63. Lower panels show the star formation history derived from the optical data, for both the differential (left, with horizontal 
dotted line indicating the past average SFR) and cumulative (right) star formation histories. The cumulative star formation history is 
calculated from the present back to 14Gyrs. Uncertainties in the lower two panels are the 68% confidence intervals, calculated from Monte 
Carlo tests including random and systematic uncertainties. Optical CMDs are restricted to the area covered by the WFC3 FOV. 
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Fig. [9I CONTINUED. — Color magnitude diagrams of the WFC3/IR (upper left) and optical (upper right) for KDG73. Lower panels show 
the star formation history derived from the optical data, for both the differential (left, with horizontal dotted line indicating the past 
average SFR) and cumulative (right) star formation histories. The cumulative star formation history is calculated from the present back 
to 14Gyrs. Uncertainties in the lower two panels are the 68% confidence intervals, calculated from Monte Carlo tests including random 
and systematic uncertainties. Optical CMDs are restricted to the area covered by the WFC3 FOV. 





10.00 1.00 0.10 0.01 Age(Gyr) 
Age (Gyr) 

Fig. Is] CONTINUED. — Color magnitude diagrams of the WFC3/IR (upper left) and optical (upper right) for KKH37. Lower panels show 
the star formation history derived from the optical data, for both the differential (left, with horizontal dotted line indicating the past 
average SFR) and cumulative (right) star formation histories. The cumulative star formation history is calculated from the present back 
to 14Gyrs. Uncertainties in the lower two panels are the 68% confidence intervals, calculated from Monte Carlo tests including random 
and systematic uncertainties. Optical CMDs are restricted to the area covered by the WFC3 FOV. 
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Fig. [9] CONTINUED. — Color magnitude diagrams of the WFC3/IR (upper left) and optical (upper right) for the target M81-DEEP within 
galaxy M81. Lower panels show the star formation history derived from the optical data, for both the differential (left, with horizontal 
dotted line indicating the past average SFR) and cumulative (right) star formation histories. The cumulative star formation history is 
calculated from the present back to 14Gyrs. Uncertainties in the lower two panels are the 68% confidence intervals, calculated from Monte 
Carlo tests including random and systematic uncertainties. Optical CMDs are restricted to the area covered by the WFC3 FOV. 
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Fig. [H CONTINUED. — Color magnitude diagrams of the WFC3/IR (upper left) and optical (upper right) for the target NGC0300-WIDE1 
within galaxy N300. Lower panels show the star formation history derived from the optical data, for both the differential (left, with 
horizontal dotted line indicating the past average SFR) and cumulative (right) star formation histories. The cumulative star formation 
history is calculated from the present back to 14Gyrs. Uncertainties in the lower two panels are the 68% confidence intervals, calculated 
from Monte Carlo tests including random and systematic uncertainties. Optical CMDs are restricted to the area covered by the WFC3 
FOV. 
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Fig. [U CONTINUED. — Color magnitude diagrams of the WFC3/IR (upper left) and optical (upper right) for the target NGC404 within 
galaxy N404. Lower panels show the star formation history derived from the optical data, for both the differential (left, with horizontal 
dotted line indicating the past average SFR) and cumulative (right) star formation histories. The cumulative star formation history is 
calculated from the present back to 14Gyrs. Uncertainties in the lower two panels are the 68% confidence intervals, calculated from Monte 
Carlo tests including random and systematic uncertainties. Optical CMDs are restricted to the area covered by the WFC3 FOV. 
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Fig. [U CONTINUED. — Color magnitude diagrams of the WFC3/IR (upper left) and optical (upper right) for the target NGC2403-DEEP 
within galaxy N2403. Lower panels show the star formation history derived from the optical data, for both the differential (left, with 
horizontal dotted line indicating the past average SFR) and cumulative (right) star formation histories. The cumulative star formation 
history is calculated from the present back to 14Gyrs. Uncertainties in the lower two panels are the 68% confidence intervals, calculated 
from Monte Carlo tests including random and systematic uncertainties. Optical CMDs are restricted to the area covered by the WFC3 
FOV. 
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Fig. [9] CONTINUED. — Color magnitude diagrams of the WFC3/IR (upper left) and optical (upper right) for the target NGC2403-HALO- 
6 within galaxy N2403. Lower panels show the star formation history derived from the optical data, for both the differential (left, with 
horizontal dotted line indicating the past average SFR) and cumulative (right) star formation histories. The cumulative star formation 
history is calculated from the present back to 14Gyrs. Uncertainties in the lower two panels are the 68% confidence intervals, calculated 
from Monte Carlo tests including random and systematic uncertainties. Optical CMDs are restricted to the area covered by the WFC3 
FOV. 
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Fig. [9] CONTINUED. — Color magnitude diagrams of the WFC3/IR (upper left) and optical (upper right) for the target SN-NGC2403-PR 
within galaxy N2403. Lower panels show the star formation history derived from the optical data, for both the differential (left, with 
horizontal dotted line indicating the past average SFR) and cumulative (right) star formation histories. The cumulative star formation 
history is calculated from the present back to 14Gyrs. Uncertainties in the lower two panels are the 68% confidence intervals, calculated 
from Monte Carlo tests including random and systematic uncertainties. Optical CMDs are restricted to the area covered by the WFC3 
FOV. 
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Fig. [U CONTINUED. — Color magnitude diagrams of the WFC3/IR (upper left) and optical (upper right) for the target NGC2976-DEEP 
within galaxy N2976. Lower panels show the star formation history derived from the optical data, for both the differential (left, with 
horizontal dotted line indicating the past average SFR) and cumulative (right) star formation histories. The cumulative star formation 
history is calculated from the present back to 14Gyrs. Uncertainties in the lower two panels are the 68% confidence intervals, calculated 
from Monte Carlo tests including random and systematic uncertainties. Optical CMDs are restricted to the area covered by the WFC3 
FOV. 
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Fig. [U CONTINUED. — Color magnitude diagrams of the WFC3/IR (upper left) and optical (upper right) for the target NGC3077- 
PHOENIX within galaxy N3077. Lower panels show the star formation history derived from the optical data, for both the differential 
(left, with horizontal dotted line indicating the past average SFR) and cumulative (right) star formation histories. The cumulative star 
formation history is calculated from the present back to 14 Gyrs. Uncertainties in the lower two panels are the 68% confidence intervals, 
calculated from Monte Carlo tests including random and systematic uncertainties. Optical CMDs are restricted to the area covered by the 
WFC3 FOV. 
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Fig. [H CONTINUED. — Color magnitude diagrams of the WFC3/IR (upper left) and optical (upper right) for the target NGC3741 within 
galaxy N3741. Lower panels show the star formation history derived from the optical data, for both the differential (left, with horizontal 
dotted line indicating the past average SFR) and cumulative (right) star formation histories. The cumulative star formation history is 
calculated from the present back to 14Gyrs. Uncertainties in the lower two panels are the 68% confidence intervals, calculated from Monte 
Carlo tests including random and systematic uncertainties. Optical CMDs are restricted to the area covered by the WFC3 FOV. 
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Fig. [H CONTINUED. — Color magnitude diagrams of the WFC3/IR (upper left) and optical (upper right) for the target NGC4163 within 
galaxy N4163. Lower panels show the star formation history derived from the optical data, for both the differential (left, with horizontal 
dotted line indicating the past average SFR) and cumulative (right) star formation histories. The cumulative star formation history is 
calculated from the present back to 14Gyrs. Uncertainties in the lower two panels are the 68% confidence intervals, calculated from Monte 
Carlo tests including random and systematic uncertainties. Optical CMDs are restricted to the area covered by the WFC3 FOV. 
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Fig. [9] CONTINUED. — Color magnitude diagrams of the WFC3/IR (upper left) and optical (upper right) for the target NGC7793-HALO- 
6 within galaxy N7793. Lower panels show the star formation history derived from the optical data, for both the differential (left, with 
horizontal dotted line indicating the past average SFR) and cumulative (right) star formation histories. The cumulative star formation 
history is calculated from the present back to 14Gyrs. Uncertainties in the lower two panels are the 68% confidence intervals, calculated 
from Monte Carlo tests including random and systematic uncertainties. Optical CMDs are restricted to the area covered by the WFC3 
FOV. 
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Fig. [U CONTINUED. — Color magnitude diagrams of the WFC3/IR (upper left) and optical (upper right) for the target SCL-DEl within 
galaxy Sc22. Lower panels show the star formation history derived from the optical data, for both the differential (left, with horizontal 
dotted line indicating the past average SFR) and cumulative (right) star formation histories. The cumulative star formation history is 
calculated from the present back to 14Gyrs. Uncertainties in the lower two panels are the 68% confidence intervals, calculated from Monte 
Carlo tests including random and systematic uncertainties. Optical CMDs are restricted to the area covered by the WFC3 FOV. 




A 
V 
00 



Redshifl (z) 
0.05 0.02 (1.01 0.005 0.002 0.001 




10.00 



1.00 0.10 

Age (Gyr) 



0.01 



o 



O 



55 0.0100 



U 



1.0000 



^ 0.1000 



0.0010 



0.0001 




10.00 



1.00 0.10 
Age (Gyr) 



0.01 



Fig. [H CONTINUED. — Color magnitude diagrams of the WFC3/IR (upper left) and optical (upper right) for the target UGC8508 within 
galaxy U8508. Lower panels show the star formation history derived from the optical data, for both the differential (left, with horizontal 
dotted line indicating the past average SFR) and cumulative (right) star formation histories. The cumulative star formation history is 
calculated from the present back to 14Gyrs. Uncertainties in the lower two panels are the 68% confidence intervals, calculated from Monte 
Carlo tests including random and systematic uncertainties. Optical CMDs are restricted to the area covered by the WFC3 FOV. 
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Fig. [9] CONTINUED. — Color magnitude diagrams of the WFC3/IR (upper left) and optical (upper right) for the target UGCA292 within 
galaxy UA292. Lower panels show the star formation history derived from the optical data, for both the differential (left, with horizontal 
dotted line indicating the past average SFR) and cumulative (right) star formation histories. The cumulative star formation history is 
calculated from the present back to 14Gyrs. Uncertainties in the lower two panels are the 68% confidence intervals, calculated from Monte 
Carlo tests including random and systematic uncertainties. Optical CMDs are restricted to the area covered by the WFC3 FOV. 
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Fig. 10. — Simulated CMDs showing the location of populations color-coded by age (upper left), initial stellar mass (lower left), and 
metallicity (upper right). The plots on the left and bottom right assume a constant star formation rate and metallicity of [Fe/H]=-1.45. 
The plot on the upper right assumes a constant star formation rate over a narrow age interval (8.9-11.2 Gyr), for a range of metallicities. 
The plot on the lower right shows the total luminosity function, along with the contributions of stars older than a given stellar age, for 
stars and ages color coded as in the CMD on the upper left. Stars younger than 2 Gyr contribute significant numbers of NIR bright 
stars. Simulations use Padova isochrones with updated AGB models from Girardi ct al. (2010), and assume photometric errors based 
on the artificial star tests for the M81-DEEP field; as we show in Melbourne ct al. (2011), the current implementation of these models 
underestimate the constribution from red core Helium burning stars compared to the data. 
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Fig. 11. — J - H vs Mh 2MASS CMD for stars in the LMC. The sequence extending up from the TRGB at Mh = -6 is dominated 
by oxygen-rich AGB stars. The horizontal sequence at Mh ^ —7.3 extending redward of J — H ^ 0.9 contains mostly carbon-rich AGB 
stars. The main sequence is the bluest vertical sequence at J — _ff ~ 0.25. The vertical sequence at J — H 0.5 is dominated by MW 
foreground stars. 
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Fig. 12. — Simulated photomet ry for a m odel galax y with a constant SFR a nd increasing metallicity (reaching Zq at the present day), 
made with the TRILEGAL code l|Girardi et al. 2005; [Girardi fc Marigo|2007[ l . The TP-AGB population is color-code d by blue /red dots 
for O-rich/C-rich stars, respecti vely, ine colors and niagmtudes ot U stars are derived from the model atmospheres bylLoidl et al.'(2001i 
(top row, panels a-c) or by Arin ger et al.| l(2009|) (bottom row, panels d-f). The colors of cool O-rich giants come fromjF'luks ct al. (19941 
model atmospheres. The left panels (a and d) show the expected H ys. J — H diagram for 2MASS filters, where C stars arc found (as 
expected) to be slightly redder than O-rich stars. In the middle panels (b and c), howev er, C stars beco me bluer than the bulk of O-rich 
stars for the WFC3/IR F160W vs. FllOW - F160W filter set,^ especially for the Aringer et al.l J2009!) models. The right panels (c and 
f) show the optical-IR color-color plots, comparable to Figure |15| For both models the the simulations include extended tails of blue and 
faint TP-AGB stars, that correspond to the ^ 15 % of TP-AGB stars caught in their low-luminosity dips along thermally-pulsing cycles, 
as well as the tails of extremely red stars at the late stages of high mass loss. 
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Fig. 13. — Comparison of CMD features in the optical (left) and the NIR (right) for IC2574-SGS. Points are color-coded according to 
their likely evolutionary phase (Blue = main sequence; Green = blue core Helium burning; Magenta = red core Helium burning; Red = red 
giant branch; Cyan = asymptotic giant branch), as deduced from the optical CMD, for stars with Mpgi4^\i[r < —2.5 that are positionally 
well-matched to stars in the NIR; unmatched or fainter stars are plotted in black. 
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Fig. 14. — Integrated Fl&OW luminosity for IC2574-SGS, integrating from the bright to the faint end, for stars classified as in Figure [13] 
Each line contains the integrated total luminosity for stars of a given evolutionary phase or younger, such that the "AGB" line includes 
the luminosity in RHeB and MS/BHeB stars as well. "Unclassified" stars are those that were not reliably matched to stars in the optical 
catalog. Note that the contribution from RHeB stars dominates the NIR luminosity for the detected stars. The vertical line indicates the 
approximate completeness limit of the data; RGB stars are likely to be the dominant contributor to the luminosity below this limit, which 
will reduce the fraction of the luminosity due to RHeB or AGB stars. Note that the "RGB" classification likely includes some contamination 
from AGB and RHeB stars as well, since such stars cannot be reliably separated within the red giant branch. 
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Fig. 15. — Optical vs infrared color for stars within IC2574-SGS, color-coded as in Figure |13| (MS=dark blue; BHeB=green; 
RHeB=magenta; RGB=red; AGB=cyan). The AGB stars span a much wider range of optical colors, due to increased variability and 
reddening from circumstellar dust. The NIR colors of the AGB fall in a tight sequence for optical colors bluer than F555W — F814VF<1.8, 
but optically redder AGB stars show increased dispersion in their NIR, colors. Only stars with magnitude uncertainties of less than 0.1 
mag are plotted. 
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Fig. 16. — (Upper Left) Color-color diagram for AGB stars from Figure |15| further split into two groups in an attempt to isolate the 
bluer, narrower sequence (blue points) from the redder population with hign dispersion in NIR, colors (red points). The NIR and optical 
CMDs of these two populations are shown on the upper right and lower right, respectively. The high dispersion selection isolates the reddest 
carbon stars (by design), and also selects against low luminosity AGB stars fainter than Afpigow = —6.5 (histograms on lower left). The 
high dispersion sub-population is also systematically bluer in FllOW — F160W at F160W magnitudes between -7 and -6.5. 
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Fig. 17. — Luminosity functions in F160W (black histogram, left axis) and median _F110Vl^ — F160W color (thick red line, right axis, 
for bins with >12 stars) for red stars ( [a] DD053; [b] DD078; [c] DD082; [d] Hoi; [e] HoIIjJf] HoII; [g] HS117; [h] 12574; [i] KDG2; ). All 
magnitudes have been extinction corrected, and distance moduli are as assumed in Table [1] The vertical dotted line indicates a fiducial 
TRGB magnitude of F160W = -5.7. 
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Fig. 1171 CONTINUED. — Luminosity functions in F160W (black histogram, left axis) and median FUOW — F16QW color (thick red line, 
right axis, for bins with >12 stars) for red stars ( [j] KDG63; [k] KDG73; [1] KKH37; [m] M81; [n] N300; [o] N404; [p] N2403; [q] N2403; [r] 
N2403; ). All magnitudes have been extinction corrected, and distance moduli are as assumed in Table [l| The vertical dotted line indicates 
a fiducial TRGB magnitude of F160W = -5.7. 
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Fig. 1171 CONTINUED. — Luminosity functions in FIQOW (black histogram, left axis) and median FUOW — F160W color (thick red line, 
right axis, for bins with >12 stars) for red stars ( [s] N2976; [t] N3077; [u] N3741; [v] N4163|_[w] N7793; [x] Sc22; [y] U8508; [z] UA292; ). 
All magnitudes have been extinction corrected, and distance moduli are as assumed in Table [l] The vertical dotted lino indicates a fiducial 
TRGB magnitude of FWOW = -5.7. 
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Fig. 18.— (Left) Absolute F160W magnitude of the TRGB as a function of FllOW - FWOW color, for data in this sample (sohd black 
circles). Also shown are loci from the isochrone models of Girardi ct al. (2008) for a 10 Gyr RGB with a range of metallicities (solid red 
line; updated to the WFC3/IR filter set), and a fit to isochrone models provided by A. Dotter (blue solid line; private communication). 
Diamonds are data from globular clusters (solid green points are from Ivanov & Borissova 2002; open yellow points are from Valenti et 
al 2004,2007). Our measurements are slightly offset from the models by roughly 0.15 magnitudes brighter in ^16014^, or 0.1 magnitudes 
bluer in color. (Right) Magnitude difference between the optical F814W (from Dalcanton et al 2009) and NIR F160W TRGB, as a 
function of the measured FllOW — FIGOW color of the TRGB. Lines and symbols are the same as in the left hand plot. Individual 
points for the Dotter isochrone TRGB values are plotted for ages of 6, 8, 10, 12, & 14 Gyr (increasing asterix sizes) and metallicities of 
[Fe/H]=-2.49, [a/Fe]=0.40, y=0.245 (magenta), [Fe/H]=-1.99, [a/Fe]=0.40, y=0.246 (blue), [Fe/H]=-1.49, [a/Fe]=0.20, y=0.246 (cyan), 
& [Fe/H]=-0.99, [o/Fe]=0.20, y=0.249 (green). 
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Fig. 19. — RGB luminosity functions for all galaxies, color-coded by their rank when sorted by the fraction of star formation in the 
most recent 1 Gyr (light lines, with redder colors indicating older mean stellar ages). Dark lines show the average luminosity function of 
galaxies, sorted by rank into 4 bins of recent star formation, with equal numbers of galaxies per bin; red, magenta, blue, and black lines 
go from lowest to highest fraction of recent star formation ((/o— iCyr) = 0.012,0.038,0.071,0.089, respectively). All luminosity functions 
are for stars found within l.Scr of the line fit to the upper 1 magnitude of the RGB, and are normalized to have the same number of stars 
within the bins between 0.25 and 0.75 magnitudes fainter than the TRGB. All magnitudes are relative to the TRGB magnitude given in 
Tabled 
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Fig. 20. — (Upper Left) FllOVK — F160W color of the TRGB, plotted against the absolute magnitude of the host galaxy in the Bt filter. 
Error bars indicate the characteristic width of the RGB near the tip (Table |3}. There is a statistically significant correlation between color 
and galaxy luminosity (Spearman rank correlation coefficient of -0.75, giving a less than 0.001% chance of occuring by chance), which is 
likely to be driven by the underlying mass-metallicity relationship for galaxies. Many of the galaxies that fall below the mean relationship 
have high recent star formation rates, which may bias the luminosity high (due to low mass-to-light ratios) or the RGB color blueward (due 
to contamination from RHeB and AGB stars, or younger mean RGB ages); ( Upper Right) Fl lOW — F160W color of the TRGB, plotted 
against the solar luminosity of the host galaxy in the Spitzer S.Gfj, bandpass l|Dale et al.||2009^ . Again, a strong color-luminosity trend is 
apparent (Spearman rank correlation coefficient of 0.82). The plot excludes ljUUYS and tic'Z'Z, which are undetected in the Spitzer 3.6/.t 
bandpass; (Lower Left) FllO W — F160W color of th e TRGB as a function the [Fe/H] metallicity inferred from the oldest bin of the star 
formation history adopted by [Melbourne et al.| ( [201 1|. There is a strong trend (Spearman rank correlation coefficient of 0.54, giving a less 
than 0.3% chance of occuring by cUance, decreasmg to 0.03% if KDG73 is excluded) driven by the fact that the color of the RGB is one of 
the strongest discriminant of metallicities in the CMD; Dashed lines indicate the fitting equations given in the text. 
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Fig. 21. — (Left) Fits to the median color of the top 2 magnitudes of the RGB, using data from Figure [T7| Colors have been randomly 
assigned to individual galaxies. All quantities have been corrected for extinction. (Right) Slopes of the fits shown in the left hand panel, 
as a function of the color of the fits at the TRGB. There are no obvious systematic trends across the wrhole sample, suggesting that the 
presence complex stellar populations can mask any intrinsic variation in the shape of the RGB with metalhcity. There are, however, 
suggestive trends at the extremes. First, the three bluest RGBs make up three of the four steepest RGBs, suggesting a possible detection 
of variations in RGB shape at very low metallicity. Second, the three fields with extremely prominent RHeB sequences (1C2574 and two 
fields in Holl) have the three shallowest RGBs, suggesting that the presence of RHeB stars are biasing the slope of the RGB by pulling the 
base to bluer colors. These same galaxies are also outliers in the luminosity- metalhcity plots in Figure [20] 



